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Abstract 



Strange particles have been a very important observable in the search for a deconfined state 
of strongly interacting matter, the quark-gluon plasma (QGP), which is expected to be formed in 
QJ ■ ultra-relativistic heavy ion collisions. We review the main experimental observations made at the 

Super Proton Synchrotron (SPS) at CERN, Geneva, and at the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory (BNL). The large amount of recently collected data 
' allows for a comprehensive study of strangeness production as a function of energy and system 

size. We review results on yields, transverse mass and rapidity spectra, as well as elliptic flow. The 
P\| I measurements are interpreted in the context of various theoretical concepts and their implications 

^ ■ are discussed. Of particular interest is the question whether strange particles are in any way 

. sensitive to a partonic phase. Finally, a compilation of experimental data is provided. 
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Table 1: Properties of the strange particles discussed in this review [TO]. Listed are the valence quark 
content, the strangeness S, the isospin, spin and parity I{J^), the mass and the decay channels which 
are most important for experiments, together with their branching ratio (B.R.) and decay length cr. 
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1 Introduction 

This review article attempts to summarize the main ideas and observations concerning strangeness in 
heavy ion reactions that have emerged during the SPS and RHIC program. Heavy ion collisions at 
these energies are believed to create energy densities that allow to explore states of matter beyond 
the deconfinement phase transition predicted by quantum chromodynamics (QCD). It separates matter 
composed of interacting hadrons from a new state of matter, the so-called quark-gluon plasma (QGP), 
where the confinement of quarks and gluons inside hadrons is removed. The order of this phase transition 
is not finally established. While for a vanishing baryonic chemical potential yUg = different lattice 
QCD studies agree on a cross over type, the situation for /Xb > is still unclear. The exact value of 
the critical temperature Tc for /ig = is also still under debate. While according to [1] Tc hes in the 
range 146 - 170 MeV, depending on the investigated order parameter, the authors of [2] find critical 
temperatures between 180 - 200 MeV. Several calculations predict that the cross over line will turn into 
a first order phase transition [31 SI thus giving rise to the presence of a critical point in the QCD 
phase diagram, while other lattice QCD investigations result in a cross over for all Hb [6]. 

The study of strange hadron production always played a special role in the investigation of QGP 
matter. Initially this was motivated by the early suggestion that an enhanced production of strange par- 
ticles, relative to p+p collisions, might provide a signature of a QGP formation. Even though strangeness 
enhancement has been established experimentally, many new theoretical developments have shed new 
light on the experimental facts. Over the years also many unexpected experimental findings have al- 
tered the viewpoint on strangeness as an observable. Strange particles have also been instrumental in 
the investigation of new phenomena that appeared in heavy ion physics (e.g. quark number scaling 
of elliptic flow). Table [1] summarizes the properties of the different strange particle species that have 
been measured in heavy ion reactions. After the bulk of the heavy ion program at the SPS has finished 
and the RHIC program has passed its 10th year, this is a good point in time to take a snapshot of the 
current situation and to summarize the main observations. 
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Figure 1: Left: The lowest order QCD diagrams for ss production ((a) qq — )■ ss, (b) gg — > ss) [7]. 
Right: Schematic picture for the reaction S + N — )■ Svr + 2K [M] . 

1.1 Strangeness enhancement as QGP signal 

The main motivation for measuring strange particles in heavy ion reactions is the expectation that their 
production rates per participating nucleon should be enhanced with respect to elementary nucleon- 
nucleon collisions, if a quark-gluon plasma is formed. Strangeness enhancement has been among the 
first signals suggested for a QGP state [HIHIE]. The argumentation by Rafelski and Miiller is based on 
two considerations. One is that in a plasma of quarks and gluons strangeness can be easily produced via 
pair production of strange-anti-strange quark-pairs. The basic processes are the fusion of two gluons or 
of two light quarks into a ss-pair (see also left panel of Fig. [1]): 

g + g s + s q + q -ir^ s + s {q = u,d). (1) 

For these reactions the Q-value corresponds to only the current mass of the produced ss-pair Qqgp = 
2ms ~ 200 MeV [10]. In contrast, the energetically cheapest way of producing strangeness in a nucleon- 
nucleon reaction is via associated production channels: 

N + N^N + A + K. (2) 

In this case the mass difference and thus the Q- value is already much higher: Qass. = itT'A + "^k — ''^n ~ 
670 MeV. As a consequence it should be much easier to generate strangeness once a plasma state has 
been formed. 

The second important point is that the equilibration times of partonic reactions, especially due 
to the gluon fusion process, are much shorter than the ones of hadronic reactions. The difference is 
especially large, if rare multi-strange (anti-)baryons are considered. In a partonic scenario with typical 
temperatures of T = 200 MeV equilibration times of Tqqp ~ 10 fm are theoretically achievable in an 
ideal gas of quarks and gluons |7]. This is on the order of the total duration of a heavy ion reaction, 
measured from the first collisions to the final freeze-out of the hadrons. However, the timespan of the 
QGP phase will be still shorter than this, so that these partonic processes might not be sufficient to drive 
the system to a complete chemical equilibrium. On the other hand, it was found in [9], that in a gas 
of free hadrons, including resonances, the typical times to reach an equilibrium state depend strongly 
on the strange particle species. While for particles with strangeness 15*1 = 1, like the kaon and A, 
chemical equilibrium might be attainable after rHG(K) ^ 30 fm, the timescales for rare (anti-)hyperons 
should be an order of magnitude longer. Following these arguments, it would thus be very difficult to 
produce multi-strange particles (S, S, Q~, in large abundances in a hadron resonance gas, while 
the presence of a QGP would be reflected in much higher production rates of these particles. In the 
latter case, the multi-strange particles would just form via quark coalescence at freeze-out and the yields 
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should be close to the chemical equilibrium expectation. A review of these theoretical considerations 
can be found in pTl [T2] . 

2 Theoretical aspects 

In the following we review the basic theoretical ideas that motivate the measurement of strange particles 
in heavy ion reactions. Beginning with the early suggestion of strangeness enhancement as a signature 
for a QGP formation, we discuss the further evolution of this idea and alternative approaches. Since 
transport and statistical models are widely employed for the interpretation of strangeness data, we 
discuss their main features in two extra sections. 

2.1 Recent developments 

Since these early calculations mentioned above, the theoretical understanding of strangeness production 
has evolved further. An overview over more recent theoretical developments is given in [13]. One 
important point is the fact that the assumption of an ideal gas, as used in the calculations in [Zl E] , is 
not valid for heavy ion reactions. As the data on elliptic flow at RHIC have shown, the matter produced 
in heavy ion collisions at high energies corresponds in its properties rather to a liquid with very low 
viscosity [TH H^l HSl HZ]. Results from lattice QCD indicate that QGP matter is characterized by the 
presence of strong correlations and bound states, which will prevail for temperatures significantly larger 
than the critical temperature [T8|, [191 EQ]. Therefore, it is inappropriate to treat partonic scattering 
processes non-perturbatively, as done in the original calculations. Attempts to take this into account, 
e.g. by hard thermal loop resummation [21], cut-off model for gluons p2], and massive, thermal gluons 
[231 SUES] 5 seem to indicate that the partonic equilibration times might rather increase [TB]. However, 
this issue is still not solved and further theoretical developments in this direction are needed. New 
insights into the problem of equilibration in a QGP might arise from approaches such as parton cascades 
[26] or the investigation of plasma instabilities [27]. Both result in a very fast equilibration on timescales 
equal or even below the expected QGP lifetimes in heavy ion reactions. Still, quantitative evaluations 
of the implications for strangeness production are missing. 

Another aspect relevant for strangeness production in heavy ion collisions is that some features of 
particle production observed at RHIC for intermediate transverse momenta pt, such as baryon- meson 
ratios and V2 scaling of elliptic flow, have raised new interest in quark coalescence as the possibly 
dominating hadronization process in this pt region. Coalescence has been suggested very early on [2S] 
and is an important ingredient of the original line of arguments to suggest strangeness enhancement as 
a signature for QGP formation [H El [9]. Recent implementations of this idea [291 EOl EU [321 [33] have 
been successful in describing the basic features of the measurements (see discussion in Section 16. ip . 

In the recent years also new mechanisms to generate strange particles in a hadron gas have been 
suggested. These include multi-meson fusion processes, such as 

Y + N^n7r + ns(Y)K, (3) 

which can additionally contribute to the yield of strange anti-baryons in a dense hadronic medium 
[31]. ns(Y) denotes the number of anti-strange quarks in anti-hyperon Y. The right panel of Fig. [T] 
shows a scheme for the reaction H + N — )■ 37r + 2K at the quark level, whose back reaction might 
give a sizable contribution to the S rates, since in a hot and dense fireball the number densities n and 
ns(Y) may be high enough. However, at RHIC energies this effect seems not to be sufficient for an 
equilibration during the lifetime of the hadronic state of the reaction [351 136] . In [37] it has been pointed 
out that the contribution from hadronic many body collisions should be strongly enhanced close to the 
phase boundary between hadron gas and QGP, due to a rapid increase of the particle densities with 
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temperature in the vicinity of Tq. The observation of equihbrated strangeness yields might therefore 
be interpreted as a direct evidence for a phase transition. Another approach that would explain a 
fast equilibration inside the lifetime of the hadronic phase is the inclusion of so-called Hagedorn states 
[38| |39] . Hagedorn states (HS) are high mass resonances produced near Tq that follow an exponential 
mass spectrum [IHl SI] and provide an efficient way of producing (anti-)baryons and kaons via reactions 
such as nvr — )■ HS — )■ nV + BB and nvr — >■ HS — )■ n'vr + KK due to their large decay widths. 

Another point of view is based on the argument that particle production via the strong interaction 
always leads to a configuration of maximum entropy, only constrained by energy, baryon number, and 
strangeness conservation [121 115 • The observed particle abundances would thus not result from 
a dynamical equilibration, but are a consequence of the hadronization process. Therefore strangeness 
yields would always be close to their equilibrium value and an enhancement in A+A collisions is rather 
due to a suppression in p+p collisions due to the much stronger effect of the conservation laws. 

Strangeness production in a hadron gas would also be enhanced, if hadron masses decrease due to 
medium modifications inside the hot and dense fireball and thus lead to reduced Q values for hadronic 
reactions. At the chiral phase transition, which lattice QCD predicts to happen before [1] or in co- 
incidence with [2] the deconfinement transition, the masses of the strange mesons should approach to 
ones of non-strange mesons, leading to a fiavor equilibrated state [IB]. In [35] it has been investigated 
within a chiral SU(3) model for the chemical and thermal equilibrium case to what extend a chiral phase 
transition is reflected in the observable particle ratios. Since this turns out to have a large effect, the 
experimentally determined particle ratios can only be understood, if the chemical freeze-out happens 
slightly below the chiral phase transition line. 

2.2 Transport models 

Hadronic transport models provide an important baseline for comparisons to the experimentally ob- 
served strangeness production. Since these models are based on hadronic degrees of freedom and string 
excitation only, they allow to explore to what extend the observed features of the data can be un- 
derstood in a purely hadronic scenario. They can therefore provide a benchmark that, to a certain 
extent, allows to judge, whether there is room for a potential partonic effect in the interpretation of 
a given experimental observation. However, when interpreting results from transport calculations one 
should be aware that these models might reach conditions (e.g. local energy densities) that cannot 
easily be reconciled with a purely hadronic scenario. But on the other hand, these type of models very 
easily allow to generate comparison data for different experimental situations, like centrality selection, 
center-of-mass energies, etc. 

Examples of hadronic transport models used in heavy ion physics for studying strangeness production 
are Ultra- relativistic Quantum Molecular Dynamics (UrQMD) [IHl HZ] , Relativistic Quantum Molecular 
Dynamics (RQMD) [18], Hadron-String Dynamics (HSD) [19]. Other models, such as A-Multi-Phase- 
Transport (AMPT) model [50] include both, hadronic and partonic, degrees of freedom. The EPOS 
model (Energy- conserving quantum mechanical multiple scattering approach, based on Partons, Off- 
shell remnants, and Splitting of parton ladders) [5T], on the other hand, follows a phenomenological 
approach that is based on a initial parton model. 

The UrQMD and the HSD model both implement binary elastic and inelastic hadronic scattering 
processes at lower center-of-mass energies and string excitation at higher ^/s^. They differ slightly in 
the number of hadronic states that are implemented. UrQMD 1.3, for instance, includes all baryonic 
resonances up to masses of 2 GeV and mesonic resonances up to 1.9 GeV, as documented by the Particle 
Data Group. Above these masses hadronic excitations are treated within a string model with meson 
formation times on the order of 1 - 2 fm/c. In the recent version 2.3 the hadron resonance spectrum has 
been extended toward a continuous distribution of resonance states for meson-baryon reactions between 
1.67 GeV < i/Scoii < 3 GeV [52]. This version includes Pythia [53j to model initial hard collisions for high 
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Figure 2: The sum of hyperons and anti-hyperons at midrapidity as a function of the number of 
participants for Pb+Pb and p+Pb colhsions at 158^4 GeV. The experimental data by the WA97 
collaboration are represented by stars, while the open symbols correspond to default UrQMD 

calculations. The full symbols denote UrQMD results derived with reduced masses or, equivalently, 
enhanced string tension [60] . 

energy reactions (a/Scou > 10 GeV). The HSD model on the other hand includes nucleons. A, N*(1440) 
and N*(1535) resonances, strange baryons and their resonances (A, E, E*, S, S*, and Q), as well as 
their anti-particles, and the 0~ and 1~ octet meson states. Above ^/s^ ~ 2.6 GeV inelastic hadronic 
collisions are described by the FRITIOF model |54j. Both model approaches are adjusted to reproduce 
the known elementary nucleon-nucleon, meson-nucleon, and meson-meson cross sections over a wide 
kinematic range as well as possible. Beyond the simple superposition of independent nucleon-nucleon 
collisions, therefore additional effects are modeled as they are present in a nuclear environment. These 
include rescattering processes between the reaction products, which can lead to strangeness exchange 
reactions like K + Y o tt + H (Y = A, S), and thus enhance the production of H. 

Generally, transport models provide a reasonable description of the production rates of particles 
with only one strange quark (kaons and A) ESI EZ] , even though there are remaining discrepancies 
in various details. However, it has frequently been shown that the observed enhancement factors of 
multi-strange particles cannot easily be understood in the context of conventional hadronic transport 
models without invoking more exotic features as, e.g., color ropes [581 ES EOl ED ES]. This is illustrated 
in Fig. [2] which shows a comparison of measured hyperons yields at 158A GeV to UrQMD results 
obtained with standard parameters. While the data for p+Pb collisions are reasonably described by the 
model, substantial differences between measurement and model are visible for central nucleus-nucleus 
reactions. In case of the Q there is an order of magnitude between data and the transport model result. 
The rates of multi-strange baryons can be matched by the model by either lowering the masses of the 
constituent quarks or increasing the string tension (k = 3 GeV/fm, compared to the standard value 
of K = 1 GeV/fm) in an artificial manner (see filled symbols in Fig. ^ [60] . So far it has not been 
possible to describe the enhancement of multi-strange baryons in a conventional hadronic scenario. 
However, multi- meson fusion processes [31], which can contribute to the yield of strange anti-baryons 
in a dense hadronic medium and might thus reduce the discrepancy between hadronic models and the 
data (see discussion in the previous section), should also be taken into account in a transport approach. 
First steps into this direction have been done within the HSD model [1^ 1^ . but a systematic and 
quantitative evaluation of this effect has not yet been performed, mainly due to the large amount of 



7 




Figure 3: Left: The number of produced partons (solid red line), mesons (long dashed green line), and 
baryon-anti-baryon pairs (short dashed blue line) as a function of time for central Pb+Pb at 158A GeV 
(impact parameter 6 = 1 fm), as calculated with the PHSD model [56l IMj . Right: The number of 
s quarks contained in anti-baryons for central Pb+Pb collisions at 158 A GeV from PHSD (solid blue 
line) and HSD (dashed red line) |56j. 

required computing resources. 

An extension of the purely hadronic transport models is the Parton-Hadron-String Dynamics (PHSD) 
transport approach, which includes also interacting partons as dynamic quasiparticles matched to re- 
produce recent lattice QCD results [561 [M]. The transition from partonic to hadronic degrees of freedom 
is described by fusion of quark-antiquark pairs or three quarks (antiquarks), obeying color neutrality, 
flavor and energy-momentum conservation. An analysis of SPS data with this model shows that up to 
40 % of the collision energy is stored in partonic degrees of freedom. The left panel of Fig. [3] illustrates 
how at top SPS energy {y/s^ = 17.3 GeV) already around 1500 partons are created in the early phase 
of the collision. This partonic phase can have a significant effect on the transverse mass distributions of 
kaons, due to repulsive partonic mean fields and initial parton scatterings, as well as on the production 
rates of multi-strange anti-baryons. As shown in the right panel of Fig. |3] the number of anti-strange 
quarks contained in anti-baryons increases by roughly a factor two in comparison to the purely hadronic 
HSD model. This approach therefore allows to estimate to which extent the presence of a partonic phase 
is needed in order to describe (multi-) strange particle production in heavy ion collisions. 
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2.3 Statistical models 

Statistical models represent another theoretical approach that plays an important role for the under- 
standing of strangeness production in heavy ion collisions. Several implementations of this model are 
frequently applied in this context, e.g. Becattini et al. [M], Braun-Munzinger et al. [66j, THERMITS 
SHARE P], and THERMINATOR The codes for latter three models are pubhcly available. 
Based on the assumption that particle yields at the end of all inelastic interactions are corresponding 
to the chemical equilibrium expectation, these type of models have generally been quite successful in 
describing measured multiplicities with a small set of parameters. The basic ingredient of these models 
are the partition functions, either grand-canonical, canonical, or micro-canonical. While the canoni- 
cal or even micro- canonical ensembles have to be used to describe small systems (e.g. p+p) the 
grand-canonical approximation provides a good description in the case of large systems such as central 
Au+Au or Pb+Pb collisions. In the grand-canonical formulation of the statistical model, the mean 
hadron multiplicities (Ni) are defined as |70] : 



The parameters are the chemical freeze-out temperature Tch, the chemical potentials fi, the volume V, 
the strange quark fugacity 7s, the spin Jj and the number of strange quarks Si of a given particle type i. 
To account for the different statistical behavior, the + sign in the denominator is valid for fermions, while 
the — sign is used for bosons. The strange quark fugacity 7s allows for the possibility that strangeness 
might not be fully equilibrated (7s < 1) and is used in some implementations of the statistical model 
[65l ITT] , while other versions assume a full equilibrium also for strange particles (7s = 1) [661 ES]- 
Other implementations allow the non-strange quarks to be out of equilibrium as well and therefore 
introduce yet an additional factor 7q [721 El] • The chemical potentials /is and /xq are determined by 
the constraint of global strangeness and charge conservation, while the baryonic chemical potential /Xb 
is a free fit parameter. Some implementations of the statistical model (e.g. [65]) directly calculate 
multiplicities according to Eq. and therefore also have to adjust the volume V. If instead particle 
ratios are fitted, as it is done in [HS], the volume parameter drops out and only Tch and /Ib remain as 
free parameter (for 7s = 1). Another aspect that has been frequently debated is whether the statistical 
model analysis should be restricted to midrapidity yields or if rather An integrated multiplicities should 
be used [651 ES]- This issue plays mainly a role for SPS energies, at which it is difficult to separate the 
central fireball and the fragmentation regions, where the physics will be different. At RHIC energies 
these regions will be well separated and for most particle species anyway only midrapidity dN/ dy values 
are experimentally available. The necessity of including 7s in the fits as a free parameter seems to some 
extend arise from the choice of An multiplicities as input. At least at SPS energies 7s comes out closer 
to one if midrapidity values are used [751 ESj . 

Figure S] shows two different fits with statistical models to heavy ion data, measured at ^/s^ = 




17.3 GeV [TB] and y/s^ = 200 GeV |77j. These examples show that the yields of all strange particles, 
including the rare Q and S"^, agree quite well with the chemical equilibrium expectation, both at the 
SPS and at RHIC. In models that include 7s as a free parameter it is found that 7s is usually larger than 
0.6 - 0.7 [78], thus indicating that strange particle production seems to be at least not far away from 
the equilibrium value. One notable exception to this behavior is the A(1520) resonance, whose yield has 
been found to be below the prediction of the statistical model by almost five standard deviations (see left 
panel of Fig. S]). Due to its short lifetime (cr = 12.6 fm) the resonance decays to a large part still inside 
the fireball. As a consequence its decay particles will be subject to elastic scattering processes in the 
medium which in turn reduce the yield of A(1520) that is reconstructible via an invariant mass analysis. 
In addition the regeneration cross section in the hadronic phase is smaller than the one for scattering 
processes of the decay particles, which results in overall suppression of the measurable resonances [7S] . 
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Figure 4: Left: The measured particle multiplicities versus the result of a statistical model fit for Pb+Pb 
collisions at y/s^ = 17.3 GeV [75]. The lower panel shows the fit residuals. Right: Measured hadron 
yields around midrapidity for Au+Au collisions at ^^/s^ = 200 GeV, together with a statistical model 
fit, resulting in the shown parameters [77] . 

This is the case for the A(1520) and the K(892) at SPS and RHIC energies [HD]- Depending on the 
suppression a lower limit for the lifetime of the hadronic phase can be derived [81]. Due to this effect 
it is suggested to exclude the resonances from statistical model fits in order to determine freeze-out 
parameters. 

It has been realized quite early that the resulting parameters T^h and /Xb lie on a single freeze-out 
curve, which follows a constant energy density over particle density of {E)/{N) ^ 1 GeV [72] (see 
left panel of Fig. |5]). This has been interpreted as a consequence of the changing composition of the 
hadron gas at chemical freeze-out, moving from a nucleon gas at low temperatures to a mainly meson 
dominated gas at high temperatures. As a result the average rest mass per particle (M) of the hadron 
gas decreases with increasing temperature, thus compensating the increasing thermal contribution to the 
energy density such that, in a non-relativistic approximation, (E)/ {N) ^ (M) +3/2 T^h ~ 1 GeV results 
[72] . Therefore, the two chemical freeze-out parameters Tch and /ia do not vary independently of each 
other, when the center-of-mass energy of the reaction is changed. If furthermore the y/s^ dependence 
of is parametrized, particle ratios at different energies can be predicted within the statistical model 
approach. In [82] the following phenomenological parametrization was suggested: 

y^B(gNN)^ , a -1.27 GeV; 6 ~ 4.3 GeV. (5) 

The resulting y/s^ dependence of Tch and is shown in the right panel of Fig. |3 Based on this 
parametrization it was shown in [82] that a maximum of the relative strangeness production is to be 
expected around y/s^ =7 — 8 GeV (see Section 11^2] for a detailed discussion). Alternative suggestions 
to parametrize the chemical freeze-out curve have been made (e.g. by requiring a fixed total (anti- 
)baryon densities of ~ 0.12fm~^ [H3], or a fixed value of the entropy density, s/T^, of approximately 
7 [SU [HS])- Other parameterizations of the energy dependence of the chemical freeze-out parameters 
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Figure 5: Left: The chemical freeze-out parameters T^h and /ig for different energies, as obtained 
from fits to data from RHIC, SPS, AGS, and SIS. The hues correspond to the freeze-out conditions 
(E) / {N) = 1 GeV and 1.1 GeV |[88]. Right: The variation of Tch and fis as a function of energy. The 
hues are the parametrization according to Eq. (|5]), using the freeze-out condition {E)/{N) = 1 GeV 
(Fig. is adapted from [S2]). 

can be found in [75]. Even though these parameterizations resuh in similar T^h - /Xb curves, they are 
motivated by different physical interpretations. For instance, fixed baryon density criterion is based on 
the picture that the freeze-out line is determined by a critical ub below which inelastic baryon-baryon 
and baryon-meson interactions cease to change the chemical composition. 

At low energies or for small systems, such as p+p, the description of particle yields with a statistical 
model has to be based on the canonical, or micro-canonical, ensemble [HI ESI EZj- In these ensembles 
strangeness production will be suppressed due to the small correlation volume in which strangeness con- 
servation has to be fulfilled [89]. Therefore, it has been suggested to model the system size dependence 
of strangeness production by a transition from a canonical ensemble to a grand-canonical one [90j . 
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3 Experimental overview 



This sections summarizes the history of the main experimental observations that have been made in 
order to estabhsh strangeness enhancement in high energy heavy ion reactions. Also, we shortly review 
the experiments at the SPS and RHIC and the experimental techniques that have been employed to 
measure strange particles. 



3.1 Basic observations 

The prediction outlined above triggered numerous experimental efforts in order to establish whether 
strangeness enhancement is realized in high energy heavy ion reactions. Early attempts at the AGS 
and SPS were limited to the measurement of particles with \S\ = 1 (i.e. kaons, A, and A), and 
to smaller reaction systems (Si+Au, 0+Au, S+S, S+Au) [9llE3[Ml[9a[95l[96l|971|98l|99]. But 
already with these data it was possible to establish that there is indeed an enhancement of strangeness 
production for kaons and A relative to pion production compared to p+p and p+A collisions ^21 [HSl 
llOOj . However, since the yields of \S\ = 1 particles measured in A+A can as well be understood in 
terms of a hadron resonance gas, including rescattering processes |55l 1101] , this observation cannot 
unambiguously be interpreted as a signature for QGP formation. Therefore, dedicated experiments 
focused in the following on the measurement of multi-strange (anti-)baryons in S+A and p+A collisions 
[TU21 nm una [ina linSl [inn Unni una inn], one of the important results of the first heavy ion period 
at the SPS was that also S~ and S"*" production is enhanced in S+W collisions with respect to p+W 
interactions |108] . With the availability of Pb beams at the SPS these studies were extended in the 
years after 1994 toward a systematic study of the system size dependence of this effect, including the 
rare fl hyperon [lll^ 1112^ 1113] . A review of the experimental results available after the first round of 
CERN experiments with Pb beams at 158 A GeV can be found in |114] . 

Figure [6] shows the strange particle enhancement in Pb+Pb collisions relative to p+Be interactions 
as observed by the NA57 experiment |112] . The strangeness enhancement factor is defined as: 



dA^(Pb+Pb) 



(iVpart) dy 



y=0 




1 dA^(p+p(Be)) 



2 dy 



(6) 



In order to account for the larger reaction volume of the nucleus-nucleus system, the measured yields 
are normalized by the averaged number of participating nucleons (iVpart)- This quantity is usually 
derived from a Glauber model calculation |115[ 1116] . A clear hierarchy of the enhancement factors for 
the different particle species is observed: While for the A (5 = — 1) the enhancement reaches values of 
3 - 4, Eg is found to be close to 10 for the E' {S = -2) and 20 for the n {S = -3) in central Pb+Pb 
collisions. Such an enhancement pattern is expected in the case of a full chemical equilibrium for a large 
system [90], which, according to the early publications [3 El E], might be interpreted as a sign for fast 
partonic equilibration in a QGP. It should be pointed out that equilibration in a partonic phase at high 
temperatures can lead to an over-saturation of strangeness production relative to the expectation of 
an equilibrated hadron gas at typical chemical freeze-out temperatures, if the transition happens very 
suddenly [117]. In any case it demonstrates that in high energy nucleus-nucleus collisions additional 
mechanisms for strange particle production are at work. The strong enhancement of multi-strange 
baryons is the most dramatic effect observed when comparing p+p to A+A collisions (other effects like 
J/ip suppression or high-pt suppression are on the order of a factor 5, while the Q enhancement reaches 
values up to 15 - 20 for central reactions). 

The NA57 results |112] are a refinement of the previous WA97 measurements |111] and thus supposed 
to supercede these data. While NA57 and WA97 use p+Be collisions as a baseline measurement, other 
experiments compare the A+A data to p+p collisions |113] I118[ 1119] . There are indications that 
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Figure 6: The enhancement of strange particle production with respect to p+Be colhsions as a function 
of the number of participating nucleons, measured by the NA57 collaboration for Pb+Pb collisions at 
158A GeV (Fig. adapted from [TT2] ). 



strangeness production is already enhanced in p+A collisions relative to p+p scaled by A^part [l^Uj 1121] 
and that a simple participant scaling is therefore violated. A projectile nucleon that collides more than 
once behaves different than one that is participating in a single binary collision (see also the discussion 
on the core corona model in Section 15. ip . Even though this effect cannot be explained by a QGP 
formation, since the created fireball is too small, and it also not sufficient to explain the enhancement 
in A+A [120j . it shows that an understanding of p+A reactions is an important prerequisite for the 
interpretation of heavy ion collisions. A consequence of this effect is that the measured enhancement 
in A+A will be less if p+A is used as a reference instead of p+p. 



3.2 Experimental techniques 

The strange particles listed in Table [1] have been measured by several experiments, both at the SPS 
and at RHIC. In Pb+Pb era at the SPS, these are the WA97, NA44, NA45, NA49, NA50, NA57, and 
NA60 experiments, while at RHIC data on strange particle production were gathered mainly by the 
STAR, PHENIX and BRAHMS collaborations (see Table [2]). A detailed summary of the available data 
measured at different center-of-mass energies can be found in section 19.11 

The reconstruction of strange particles is done via two basic experimental methods, which are 
employed by the experiments listed above. On one side there is the direct reconstruction as a charged 
particle track (K^) and the indirect reconstruction via their weak decay topology (Kg, hyperons, and 

via kinks). In the following, a short description of the basic methods is given and their strengths 
and associated difficulties (e.g. required corrections) are discussed. 
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Figure 7: The separation of kaons to pions and protons as achievable by specific energy loss (dE/dx) 
measurements for a resolution of 4% as obtained by the NA49 experiment [123] . 



3.2.1 Charged kaons 

In the case of charged kaons, a direct identification of the particles can be performed. Their trajectory 
can be reconstructed and their identity derived from either their specific energy loss (dE/dx) in the 
detector material, their time-of- flight (TOF), the kink decay topology, or their Cerenkov signal. 

The main advantage of using the specific energy loss for particle identification is that it is usually 
measured with the same device that is used for tracking (e.g. time projection chambers (TPC) or silicon 
detectors) and that therefore the same acceptance can be covered. This allows to identify particles in 
large regions of phase space. The charged particles traveling through the active volume of the detectors 
lose energy by atomic collisions in the gas (TPC, drift chamber) or silicon (Si-detectors). The released 
charge depends only on the velocity of the particle, as described by the Bethe-Bloch formula 



dx A/32 




P ] ^ 2 



(7) 



and thus allows to separate different particle species experimentally, if their momentum is known. K 
is a constant, z the charge of the incident particle, m^. the electron mass, Z the atomic number and 
A the atomic weight of the detector medium, /3 the velocity of the incident particle in units of c, 
Traax the maximum kinetic energy which can be imparted to a free electron in a single collision, I the 
mean excitation energy of the detector medium in eV, 7 = (1 — I3'^)~^/'^^ and 5(/37) the density effect 
correction [1^. At lower momenta (p < 3 — 4 GeV for pions), where the specific energy loss follows 
the l//3^ dependence of the Bethe-Bloch curve, the kaons can be easily separated from protons and 
pions. Therefore the dE/dx method is very well suited for kaon identification in collider experiments 
up to momenta of p < 0.7 GeV/c |122] . Above this momentum the kaon band crosses the pion band 
and a clear separation is difficult. In fixed target experiments at higher energies, where due to the 
additional Lorentz-boost particles have to be measured at higher laboratory momenta, the dE/dx 
method becomes more involved. Because in the high momentum region the specific energy loss is 
described by the relativistic rise region of the Bethe-Bloch curve, where the bands of pions, kaons, and 
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Table 2: Overview on experiments measuring strange particles at the SPS and at RHIC. 



Experiment 


Particles 




Reaction systems 


v^(GeV) 


4:71 


NA44 






Pb+Pb 


17.3 





NA57 (WA97) 


K°, A, s, n 




p+A, Pb+Pb 


8.7, 17.3 





NA45 


e+ + 




Pb+Pb 


17.3 





NA49 


K±, 0, A, E, 


^] 


p+p, C+C, Si+Si, Pb+Pb 


6.3, 7.6, 8.7, 12.3, 17.3 


yes 


NA50 


yU"^ + yU" 




Pb+Pb 


17.3 


— 


NA60 


— i- /i+ + /i" 




In+In 


17.3 


— 


STAR 


K±, 0, A, S, 




p+p, d+Au, Cu+Cu, Au+Au 


9.2, 62.4, 130, 200 




PHENIX 


K±, 0, A 




Au+Au 


130, 200 




BRAHMS 






p+p, Au+Au 


62.4, 200 


yes 



protons are relatively close together. But also in these cases the dE/dx signal has been successfully 
employed for kaon identification. For instance the NA49 experiment was specifically designed to reach 
a dE/dx resolution of 4% |123j (the resolution for STAR is better than 8% |T24]), which is sufficient 
to de-convolute the measured dE/dx spectra and to extract charged kaon yields on a statistical basis 
|125[ 1126] . Figure [7] shows the separation power that can be achieved with this resolution. 

The time-of-fiight method, in many cases combined with a dE/dx measurement, can provide a clean 
identification of charged kaons in the momentum range p < 2.0 GeV/c [1251 11261 11271 11281 1129] . To 
apply this method, the path length and the momentum of a given particle has to be known, as well 
as the precise time-of-fiight. The required time resolution typically has to be below 100 ps, which can 
easily be achieved with scintillator based TOF detectors in combination with a fast start detector. Since 
the TOF method requires a separate detector, the acceptance region of the TOF measurement is in 
most cases, due to cost reasons, smaller than the one available for tracking. Modern techniques, such as 
resistive plate chambers (RPC), allow to build TOF detectors with a large acceptance. These detectors 
can cover, e.g., the full 27r azimuth of the TPC in the STAR and ALICE experiments. 

Sometimes an alternative way of reconstructing charged kaons is followed, which exploits the topol- 
ogy of a weak decay, namely K"*^ fi^ + or n^ + 7r°. Since the neutral particles are not 
detected, the reconstruction algorithms try to identify this decay by looking for kinks in the tracks of 
charged particles. This method is complementary to dE/dx or TOF measurement and can provide an 
independent way of particle identification with different systematic effects |130] . However, there are 
substantial sources of background, like charged hyperon and pion decays, as well as multiple scattering 
in the detector material jlI4j . The method is limited to particles that decay inside the sensitive area 
of the tracking devices which reduces the available statistics. 

3.2.2 Hyperons and Kg 

The other method, which is used for the majority of strange particles, is based on the analysis of their 
decay topology. This method requires the simultaneous reconstruction of 2 (Kg, A, A) or 3 (S~, S+, 
VL^) charged tracks. Figure [8] illustrates the topologies of a A, respectively decay. In order 
to search for V° candidates (i.e. Kg, A, and A), all oppositely charged tracks are combined into pairs 
and the distance of closest approach (DCA) between them is searched for. The position of the DCA 
gives the possible secondary vertex position. Since weak decays are characterized by long lifetimes of 
the mother particles, a large fraction of their decay vertices has a substantial spatial distance to the 
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30 cm 



Figure 8: Left: The topology of the decay A — > p + vr^. The sohd hnes correspond to the reconstructed 
particle tracks of the proton and the negative pion, while the dotted ones represent their extrapolations 
through the magnetic fields toward the target plane. Shown is the distance of closest approach (DCA) 
between the two charged tracks. Right: Illustration of the reconstruction of the decay Q~ — )■ A + K~ 
within the NA57 setup (not to scale) |112j . 

interaction vertex. Therefore, the requirement of a minimal distance between both allows to reduce 
the combinatorial background, which to a large extend results from combining primary tracks. As 
the corresponding decay lengths are on the order of centimeters (see Table [1]), the secondary vertex 
resolution has to be in the same order or better for collider experiments. The requirements are less 
stringent in the case of fixed target experiments, since the average decay lengths are larger due to the 
additional Lorentz boost. In many cases (e.g. E895 [m], NA49 [58], STAR [1321 [133], and PHENIX 
|134j ) TPCs or drift chambers provide sufficient spatial resolution for this purpose. Other experiments 
(e.g. NA57 [112] ) use silicon detectors as a tracking device, which provide much better vertex resolution 
in the sub-millimeter region and a higher rate capability at the expense of a reduced acceptance. In large 
collider experiments both techniques, TPC and silicon, are combined (e.g. ALICE |135] and STAR). 
By applying additional topological cuts (e.g. on the pointing angle of the mother momentum vector 
relative to the main vertex position) or using additional particle identification on the daughter tracks, 
the combinatorial background can be further reduced. For the particle pair assigned to a given V° 
candidate the invariant mass is calculated based on a mass hypothesis by assigning either the pion mass 
to both tracks (Kg) or the pion and the proton mass (A) to the positively and negatively charged track, 
respectively. From the resulting invariant mass distributions the yields of the corresponding mothers 
particles can be extracted. Typically a mass resolution around 4 MeV is achieved. 

A and Kg can also be reconstructed without measuring the decay vertex position, as has been 
demonstrated by PHENIX |134] and NA45 [136] . The drawbacks of this approach are an increased 
combinatorial background and a poorer mass resolution, while the benefit is a somewhat reduced sen- 
sitivity to the accuracy of the detector simulation. 

A (A) candidates from this first step can then be combined with a third charged track to form S~ 
(S"*") and (^^) candidates. Since the A are in these cases daughter particles, it is usually required 
that their momentum vector does not point back to the interaction vertex, in contrast to the case of the 
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primary A. This method has been successfully applied to heavy ion reactions to measure systematically 
the production of the rare Q~ and [591 1112[ 1133] , although usually the sum Q~ + is shown in 
order to increase the statistical significance. 

3.2.3 Corrections 

The different detection techniques for strange particles differ in their necessary corrections. The ones 
for the direct measurement of charged kaons are mainly determined by the single track acceptance 
and reconstruction efficiency, which is usually quite high (typically > 90 % at higher momenta), thus 
resulting in relatively small correction factors. In contrast, the reconstruction of the decay topology 
requires the measurement of two charged tracks. Therefore any acceptance losses and inefficiencies enter 
quadratically. Additionally, the algorithms for the secondary vertex reconstruction usually cause losses 
which have to be corrected, mainly due to the cuts applied to reduce the combinatoric background. As 
a consequence, the corrections for this method are higher than for the single track measurements and 
have to be carefully determined by a Monte Carlo simulation (e.g. for A the efficiencies vary typically 
between 5 % and 50 %, depending on momentum). On the other side, the particle identification of the 
topological method is clean and unambiguous. 

Another correction that is necessary for some particles is the correction for feed-down from weak 
decays. This is important for A and A, where part of the measured yield is originating from the decays 
S~(H"'") — )■ A(A) + 7r~(7r"'') 0. Due to the long decay lengths of the weak decays, the fraction of accepted 
feed-down particles strongly depends on the detector geometry, the reconstruction algorithm, and the 
applied cuts (typical values for the size of the feed-down correction are 5 % to 20 % for A, for A it 
can be even higher). Therefore, this correction has to be applied at the analysis stage and cannot be 
calculated properly afterward without exact knowledge of the experimental conditions. This argument 
does not apply to strong and electro-magnetic decays, which all happen close to the reaction zone. 
There, all decay particles will be included in the measured yields and can thus easily be subtracted, if 
necessary. For instance, the measured A multiplicities usually refer to the sum A + since A from 
the decay S'^ — A + 7 cannot be separated experimentally from the primary ones without performing 
an additional photon measurement. 



^The weak decays Q {^'^) — >■ A(A) + K (K"'") are usually ignored, since their contribution is marginal. 
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Figure 9: The p/p, A/A, , and f2+/i7 ratios around midrapidity in central Pb+Pb and Au+Au 

collisions as a function of ^/s^. Shown are data by the SPS experiments NA49 EH 1137] . NA44 
[188] and NA57 [1121 [I39], from AGS ED Ea ES] and RHIC experiments [ISS UMl 



4 Energy dependence 

The large amount of data at different center-of-mass energies that have been accumulated at the SPS and 
at RHIC over the recent years allows, together with AGS measurements from low energies, to compile 
a relatively complete picture of the energy dependence for most strangeness related observables. Of 
special interest is the ^/s^ evolution of the observed strangeness enhancement. Also, several intriguing 
structures in the energy dependences have emerged, for instance in the K^/tt^ ratio. The (p nieson, 
which might be directly sensitive to the partonic phase of the collision, deserves a special attention and 
is treated in a separate subsection. We also discuss the ^/s^ dependence of transverse and longitudinal 
spectra. 

4.1 Strangeness enhancement at different energies 

In order to answer the question whether the observed strangeness enhancement can really be caused by 
a fast equilibration in a partonic phase, the investigation of its energy dependence provides important 
information. Following this argument one should observe an onset of strangeness enhancement when 
going from low center-of-mass energies toward higher ones. At lower energies, where the energy densities 
should be too low to form a deconfined state, any enhancement relative to elementary collisions can 
only be due to secondary hadronic interactions, which are, at least in the case of multi-strange baryons, 
unlikely to drive the system toward chemical equilibrium. It has been pointed out quite early that the 
enhancement for kaons is, in contrast to naive expectations, rising toward lower energies |15H 1152] . A 
similar observation has been made for A [58] , while for multi-strange particles the current experimental 
situation still is unclear. Data on multi-strange particles for heavy ion collisions at lower energies are 
scarce and have large errors, while p-f-p reference data are essentially not existent. While first attempts 
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Figure 10: The enhancement factor for midrapidity yields measured in central Pb+Pb (Au+Au) 
collisions at y/s^ = 17.3 GeV (200 GeV) versus the number of strange valence quarks for mesons 
(left) and baryons (right) [118] . Shown are data by the NA49 collaboration (filled circles) [5S | 15^ 1 ri25[ 
HiTl Hm [US], by the NA57 collaboration (open circles) [112], and by the STAR collaboration (filled 
stars) |119[ I122[ I150j . NA49 and STAR use p+p collisions as baseline, while in the case of NA57 Es is 
calculated relative to p+Be collisions. 

into this direction have been done for heavy ion collisions |13H 1153] . there is still a lack of reference 
data from p+p or p+A collisions at the same energies. Threshold effects play an important role at 
very low energies and complicate any interpretation. Nevertheless, it is interesting to observe that S~ 
measurements at very low energies are close (Au+Au at 6 A GeV, E895 collaboration [131] ). or even 
above (Ar+KCl at 1.76 A GeV, HADES collaboration |153] ). the expectation for an equilibrated hadron 
gas [821 [TM]. 

At higher energies, a clear trend in the energy dependence of strangeness enhancement emerges. 
The comparison of the hyperon enhancement factors measured at top SPS energy (y^s^ = 17.3 GeV) 
and top RHIC energy {y/'s^= 200 GeV) shows that the enhancement for baryons is actually higher 
at SPS. As shown in the right panel of Fig. [TOl the enhancement factor Es (Eq. ([6])) is larger for S~ 
and at y/s^ = 17.3 GeV than at y/s^ = 200 GeV. However, for A and mesons (K^, K^, 0) the 
enhancement is essentially the same at both energies (left panel of Fig. [TOj) . Preliminary data by the 
NA57 collaboration indicate that in central Pb+Pb collisions at y/s^ =8.7 GeV the S~ enhancement 
is again larger as compared to p+Be collisions than at y/s^ = 17.3 GeV [155] . The main reason for 
this behavior is that the strangeness production in p+p (p+Be) collisions is rising faster with energy 
than in nucleus-nucleus reactions. For instance, the S+Ztt"*" ratio changes by a factor of ~ 3.8 in p+p 
and of only ~ 2.9 in central A+A between ^s^^ = 17.3 GeV and 200 GeV [SH |122l [I33l ITi71[T56] . 
This is even more pronounced for the /Tr~ ratio, since for central nucleus-nucleus collisions this ratio 
is decreasing between the two energies by roughly 20% (see Fig. [T^ . while at the same time it is rising 
in p+p collisions by 25%. In the context of statistical models the suppression of strangeness production 
is understood as a consequence of the small available volume, requiring a canonical ensemble, which 
results in a strong reduction of the available phase space for strange particles. Along with the center- 
of-mass energy the overall multiplicity in p+p reaction increases, and thus the effective volume, with 
the result of a less effective phase space suppression. Proton-proton collisions therefore slowly approach 
the chemical equilibrium values for large systems (which then can be described by a grand-canonical 
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Figure 11: Left: The energy dependence of the Wroblewski factor As as determined via fits of a 
statistical hadron gas model to multiplicities measured in elementary (e+e~, K^p, 7r+p, pp, pp) as well 
as central heavy ion collisions [157] . The dashed line is a prediction for a fully equilibrated hadron 
gas, the solid line an interpolation of fits including an additional strangeness suppression factor 7s 
(Fig. is adapted from |157] ). Right: The energy dependence of the strangeness undersaturation factor 
7s at chemical freeze-out [78]. The dashed line is a phenomenological parametrization in the form 
7s = 1 — 0.606 exp |— 0.0209 ^ A |, where A is the atomic mass number, fulfilling 7s — )■ 1 for 
^s^^oo [78]. 



ensemble), such as heavy ion collisions. 

The A deviates from the general trend, because its enhancement significantly increases with center- 
of-mass energy. But in this case the strong change of the net-baryon density in A+A collisions with 
energy and the consequently strongly changing anti-baryon/baryon ratio is dominating the energy 
evolution of the A enhancement. The sensitivity of a given anti-baryon to this effect depends on 
its strangeness content. As a consequence a clear hierarchy of the anti-baryon/baryon ratio R is visible, 
Ripf^ /Vt~) > R(E'^ /E^) > i?(A/A) > i?(p/p), and the energy dependence of the ratios gets significantly 
weaker with increasing strangeness content (see right panel of Fig. [9]). 

Another way of looking at the energy dependence of strangeness enhancement is studying the Wrob- 
lewski factor As [65] . It is defined as the ratio of the mean multiplicities of the newly produced valence 
quark-antiquark pairs [158] : 

- (uu) + (dd) 

Since usually not all hadron multiplicities for a given reaction system are measured, the determination of 
As has to rely on some input from models. In Fig. [11] As has been extracted from fits with a statistical 
hadron gas model for elementary collisions, as well as for central heavy ion reactions [157] . The As 
values determined for nucleus-nucleus collisions are found to be a factor 2 higher than for elementary 
collisions, even at the lowest SPS energy {y/s^ = 6.3 GeV). There is an indication for a rise of As at 
lower energies (y^s^ < 6 GeV), but from then on As stays at a constant value of ~ 0.5, well above the 
value of ~ 0.25 found for elementary collisions. This demonstrates again that strangeness production is 
globally enhanced in heavy ion reactions, already at very low energies. However, it remains to be seen 
whether this is not only true for the global strangeness production, but as well for the rare multi-strange 
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Figure 12: Left: Lines of constant Wroblewski factor As in the Tch — /ie plane. The thick sohd hne 
without labehng represents the chemical freeze-out curve defined by {E)/{N) ~ 1 GeV |82]. Right: 
The contributions to the Wroblewski factor from strange baryons (dotted line), strange mesons (dashed 
line), and with hidden strangeness (dash-dotted line) in the full equilibrium case (i.e. 7s = 1) |52]- The 
sum of all contributions is given by the solid line. (Figs, are adapted from |H2]-) 



baryons. 

In the statistical model approach by Becattini et al. [78] the strangeness undersaturation factor 
7s is a free parameter in the fits to the particle yields. The right panel of Fig. [11] shows the energy 
dependence of 7s for nucleus-nucleus collisions between AGS and RHIC energies. A slow increase with 
y/s^ from 0.7 - 0.8 at SPS toward a value of 7s = 1 at RHIC is observed. This would indicate that 
strangeness production is not yet fully equilibrated at the SPS, even though it is already very close to 
the full equilibrium case. 

Under the full equilibrium condition (i.e. 7s = 1) the statistical model results in a well defined 
connection between a given chemical freeze-out position in the T^h — Hb plane and the Wroblewski 
factor. The left panel of Fig. [T2] shows the lines of constant As in the Tch — /^b plane [S2] • Based on 
the parametrization of the y/s^ dependence of /Zb, as e.g. given in Eq. ([5]), and the (E) / (N) ^ 1 GeV 
freeze-out criterion (solid line in the left panel of Fig. [TT]) . which connects Tch to a given /ig, a unique 
curve for the energy dependence of As can thus be constructed [82]. This curve is shown as solid line in 
the right panel of Fig. Upland is characterized by a distinct maximum of relative strangeness production 
around ^/s^ of 6 - 8 GeV (corresponding to /i-e ~ 500 MeV and Tch ~ 130 MeV). The dotted line in the 



21 



left panel of Fig. [TT] is based on a similar approach, but uses different parameterizations of the energy 
dependence of the chemical freeze-out parameters [157] . It follows from these considerations that the 
relative contributions to As change with energy. While at lower energies, where /Ib is large, the produced 
strange quarks are to a large extent contained in strange baryons (dotted line), the relative strangeness 
production at higher energies {y/s^ > 6 GeV, corresponding to /ig < 500 MeV) is dominated by strange 
mesons (dashed line). The maximum in As at lower energies is thus mainly due to the contribution of 
the strange baryons. Also strange mesons contribution exhibits a maximum, which is, however, less 
pronounced and shifted toward higher energies. This expected maximum of relative total strangeness 
production is reflected in the energy dependence of various strange to non-strange particle ratios, as 
discussed in the following sections. 
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Figure 13: The energy dependence of the (K"'")/(7r"'") ratio (left panel) and the (K~)/(7r~) ratio (right 
panel) measured in central Pb+Pb and Au+Au collisions (solid symbols) |125|. I126i I159[ I160[ I161[ I162[ 
11631 1164^ I165[ I166j compared to results from p+p(p) reactions (open symbols) [ 15H I167j and from the 
transport models HSD (solid line) ^ |S3], UrQMD (dotted line) ^46^, and RQMD (long dashed hue) 
pS] and from a statistical hadron gas model (short dashed line) [721 ES] • 



4.2 Structures in the energy dependence of particle yields 

Apart from the search for an onset of strangeness enhancement, the measurement of strange particle 
abundances in heavy ion collisions at different center-of-mass energies revealed new insights into the 
physics of strangeness production. Data from several experiments at the AGS, SPS, and RHIC allow 
to compile a comprehensive overview for central A+A collisions. In the following, only some of the 
most interesting features seen in the energy dependences shall be discussed. A compilation of the ^/s^ 
dependences of the midrapidity yields of all strange particles is given in section 19.11 (see Figs. HQ] and Hi] 
and Tables IMID. 

4.2.1 The K/tt ratios 

One of the most prominent observations in this context is a pronounced maximum in the energy de- 
pendence of the K^/Tr"*" ratio around y^s^ = 7 GeV |126j . as shown in the left panel of Fig. [13] The 
K~ /tt^ ratio, in contrast, exhibits a rather smooth energy dependence with only a slight indication for a 
kink at the position of the maximum in the /7t~^ ratio (right panel of Fig. [T3|l . Transport models like 
RQMD [3S], UrQMD ^B], and HSD |1 ^ [55 [[TBS] are generally not able to describe the structures seen in 
the data. For instance, they just predict localized turnover in the energy dependence of the K^/tt"'" ra- 
tio, but not the observed maximum. Similar discrepancies are present for the K~ /tt~ ratio, although 
they are not so visible in this case, due to the less pronounced structure in its energy dependence. 

On the other hand, a structure like the observed one has been predicted in |169] for the energy 
dependence of the strangeness to entropy ratio as a consequence of the onset of deconfinement. This 
prediction is based on the statistical model of the early stage [ 170] , which is an extension of Fermi's 
statistical model of particle production that takes into account two different states of matter, a confined 
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Figure 14: Left: The energy dependence of the rapidity densities dN/dy around midrapidity of K"*" and 
A, divided by the ones of tt"*", respectively vr", for central Pb+Pb and Au+Au collisions. The solid 
line is the result of a statistical model calculation. The dotted line gives the K'^/tt'^ ratio including 
the additional effect of higher mass resonances. The dashed lines show the energy dependence of the 
temperature Tch (upper panel) and of the baryonic chemical potential /Xb (lower panel) [77]. Right: 
The energy dependence of the rapidity densities dN/dy around midrapidity of A, A, S~, and S"*" 
divided by the total pion rapidity densities (tt = 1.5 (tt" + 7r+)) for central Pb+Pb and Au+Au collisions 
[SSl[II3[I3ni[ffl[I3Sl[IMl[I3HlM[ffl[IlSl[^ Please note that the vertical error 

bars correspond to the statistical errors only. Also shown are results for the transport models UrQMD 
(dashed line) |16] and HSD (dash-dotted line) jl9l |63] , as well as a statistical hadron gas model (solid 
line) [75]. 

(hadron gas) and a deconfined (QGP) state. The initial state is described by Fermi-Landau conditions, 
which allows to define the energy density at a given center-of-mass energy. At a critical temperature Tq a 
first order phase transition is assumed. This temperature corresponds via the energy density to a certain 
center-of-mass energy and is adjusted such that the transition occurs around y/s^ ^7 — 8 GeV. The 
model assumes full statistical equilibrium in the initial stage of the reaction. Strangeness and entropy 
are conserved throughout the evolution of the system. Once the transition energy is passed, a change in 
the energy dependence of the strangeness to entropy ratio is expected due to the change of mass of the 
strangeness carriers and of the available number of degrees of freedom. At low energies strangeness is 
confined in hadrons with high mass, leading to a steep energy dependence of the strangeness to entropy 
ratio. At high energies, above Tq, the mass of the strangeness carriers is assumed to be reduced to the 
mass of the strange quarks, resulting in an almost flat energy dependence of this ratio. The charged pion 
multiplicity is assumed to be proportional to the entropy production. Since the K"*" are carrying almost 
half of the produced anti-strange quark^ and are thus directly sensitive to the total (anti-)strangeness 
production, the /tt~^ ratio should exhibit the same behavior as the total strangeness to entropy ratio. 

^The other half is mainly contained in K*^, due to isospin symmetry, while strange anti-baryons do not contribute 
significantly at lower energies. 
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Figure 15: Left: The energy dependence of the {Q)/{tt) ratio in central Pb+Pb colhsions [22] • Also 
shown are results for the transport models UrQMD (dashed line) [H] and RQMD (hashed box) |172] . 
as well as of statistical hadron gas models, one with strangeness undersaturation factor 7s (dotted line 
with empty boxes) [76] and one without (solid line) [82]. Right: the energy dependence of the fl~ and 
Cl^ rapidity densities dN/dy around midrapidity for central Pb+Pb |139] and Au+Au |132] collisions. 

Even though this interpretation is highly debated and relies on simplified assumptions (e.g. ideal gas 
equations), the structure in the K^/7r+ ratio might nevertheless be indicative for a sudden change in 
the nature of the reaction system produced in this energy region. 

The statistical hadron gas approach, using a parametrized y/s^ dependence of the chemical freeze- 
out temperature Tch and the baryonic chemical potential /Xb (see Eq. (^), as based on fits to the 
measured particle ratios [721 [82], also results in a maximum of the /7r~^ ratio at approximately the 
right position, although its width is clearly wider than the data (see short dashed curve in left panel 
of Fig. [T3|) . However, a recent extension of this model that includes the scalar a meson and high 
mass resonances (m > 2 GeV) [77], following the early ideas of Hagedorn [IDl HI], leads to an improved 
description of the data (left panel of Fig. [T^ . In [77] it has therefore been argued that the special form of 
the energy dependence of the /7t~^ ratio would be a consequence of reaching the limiting temperature 
predicted by Hagedorn |173j and therefore indirectly indicate the presence of a quark-gluon plasma 
phase. 



4.2.2 Other particle ratios 

Similar maxima as in the case of the K'^ /7t~^ ratio, although not as sharp, can be observed in the ratios 
of strange baryons to pions (panel (a) and (c) of Fig. UM [58J. While the A/vr ratio is relatively well 
described by transport models, the S~/7r ratio is clearly under-predicted by the HSD model [IHl 1^ . 
Generally, hadronic transport models are not able to match the measured abundances of multi-strange 
(anti-)baryons (see Fig. [MT d) and left panel of Fig. [T5|) . The primary production mechanism of string 
fragmentation does not favor the generation of particles with multiple strangeness without going to 
unrealistic parameter settings (see discussion in Section [2^2]) and subsequent rescattering processes add 
only little to the final multiplicities. Multi-meson fusion processes might help to approach the data, 
but the currently available transport codes do not yet include this feature. 

Statistical hadron gas models generally provide a relatively good description of multi-strange par- 
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Figure 16: The K~ /K~^ ratio as a function of the p/p ratio [1 74] . The data are taken from several 
experiments \l'25\ 1126^ 1137^ 1140^ 1166^ 1174^ 1175] and are measured for central collisions at different y/s^ 
and rapidities. Also shown are various model predictions [IHl SZl EOl 1176] . 



ticle^ [ZOl [75], including the f2 (left panel of Fig. [15]). So far, all measurements of strange particle 
multiplicities in heavy ion reactions are in reasonable agreement with the expectation for an equi- 
librated hadron gas at all energies, although in some implementations of this model, which use Air 
integrated multiphcities as input, the additional strangeness undersaturation parameter 7s is needed to 
provide a proper fit to the data everywhere [76] (see discussion in Section |23]) . There is no indication 
yet for any significant departure from the chemical equilibrium values at lower energies |131[ 1153] . 

In the statistical model approach the maxima of the A/vr, H~/7r, and K^/tt ratios are due to the 
specific dependence of the baryonic chemical potential /Ib on the center-of-mass energy. This results in 
a maximum of relative strangeness production around y/s^ = 7 GeV [82j. The right panel of Fig. [T2] 
shows how, according to this model, the Wroblewski factor As changes with energy. Especially the 
contribution from strange baryons exhibits a distinct maximum around y/s^ = 4 GeV, which in the 
data is most clearly visible in the A/tt ratio (Fig. [T^ a)). The position of the maximum is expected 
to shift toward higher energies with increasing strangeness content of the baryon, and the maximum 
should be less pronounced for the S~ and the Q~ [177] . The observed peak in the H^/tt ratio seems to 
be smaller in the NA49 measurement [SB], while the NA57 data suggest a stronger energy dependence 
|139] (see Fig. [TW c)). For the Q there is no evidence for a maximum, but the current data are limited 
in significance and coverage (see Fig. [15]) . 

Figure [TBI shows the K~/K'^ ratio as a function of the p/p ratio [174] . The data points measured at 
RHIC energies by the BRAHMS collaboration {y/s^ = 62.4, 130, and 200 GeV) have been measured 
at different rapidities, while the SPS and AGS data correspond to midrapidity. Generally, a universal 
dependence of the K~ /K~^ on the p/p ratio is observed. Also included are predictions from two different 
transport models (UrQMD [IHl HZ] and AMPT \5U\ 1176] ). Even though they roughly follow the trend, 
they do not result in the universality observed in the data. In the statistical model approach, on the 

■^Please note that concerning the hyperon muhiphcities there are stiU remaining discrepancies between the experiments 
NA49 and NA57 at the SPS (Figs. [TH and HI]), which exceed the quoted systematic errors. These are present for all particles 
at 158A GeV and for A and S~ at 40A GeV. Possible reasons might be difficulties in correcting trigger inefficiencies 
(NA57) or an insufficient correction of the multiplicity dependent reconstruction efficiency (NA49). But despite extensive 
discussion between the collaborations the origin of these discrepancies were never finally resolved. 
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Figure 17: Left: The (p/n ratio at midrapidity in central nucleus-nucleus collisions as a function of the 
center-of-mass energy [149J. Also shown are the predictions of a statistical hadron gas model [75j and of 
the transport model UrQMDl.3 [46J. Right: The total multiplicity in central nucleus- nucleus collisions 
as a function of the center-of-mass energy [149] . The filled circles, connected by a solid line, represent 
results from a fit with a statistical hadron gas model that includes a strangeness undersaturation factor 
7s [ZH]- The dotted curve shows the prediction of the UrQMD model [16]. (The Figs, are adapted from 

m-) 



other side, the p/p ratio is defined by the baryo-chemical potential /ig and the chemical freeze-out 
temperature Tch as |178] 

£ = exp (^) . (9) 

Thus, the p/p ratio, measured in the given rapidity range, should provide a good estimate of the local 
/is- Under the assumption that the freeze-out temperature T^h is constant, the p/p ratio would therefore 
fully determine the chemical composition of the system. A fireball observed at high rapidities and high 
energies (e.g. ^/s^ = 200 GeV) would be equivalent to one at midrapidity and lower energies (e.g. 
y/s^ =17.3 GeV), and the measured particle ratios should be the same, as shown for the K'/K"*" ratio 
in Fig. [ini A detailed analysis |179] shows that both parameters, /Zb and Tch, depend on rapidity, but 
can be related via the freeze-out condition {E) / (N) ^ 1 GeV [72] . However, the dependence of /Ib on 
rapidity is found to be stronger at SPS energies than at RHIC. 



4.2.3 The Meson 

Among the strange particles, the (j) meson has been of special interest. Due to its ss valence quark 
structure it is strangeness neutral and should therefore not be sensitive to any hadro-chemical effects, 
as they are described in statistical hadron gas models by a strangeness undersaturation factor 7s or 
the canonical suppression for small systems. However, if the production was dominated by quark 
coalescence from a partonic phase of the reaction, its sensitivity would be even more pronounced. In 
this case, any strangeness undersaturation should affect the rather as 7^. On the other side, in a 
hadronic phase mesons can be produced by kaon coalescence, K"*" -|- K~ — )■ 0. 

The energy dependence of the 0/7r ratio, measured around midrapidity in central nucleus-nucleus 
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Figure 18: Upper left: The energy dependence of the (f>/K~ ratio in nucleus-nucleus (full symbols), 
e"'"+e~ (open triangles), and p+p (open circles and star) collisions |180j . Lower left: The system size 
dependence of the 0/K~ ratio in different collisions systems. The dashed line shows scaled results for 
the transport model UrQMD jlB]. Right: The energy dependence of the (p yield, normalized by the 
product of the and K~ yield (top), the square of the yield (middle), and the product of the 
yield with the sum of the K~ and A yields (bottom). All yields are measured around midrapidity 
for central Pb+Pb or Au+Au collisions. Shown are only statistical errors. The data are taken from 

collisions, is summarized in the left panel of Fig. [T7] [HSj. This ratio increases monotonically with 
center-of-mass energy without any significant structure. The data are compared to a statistical hadron 
gas model, assuming full strangeness equilibration (i.e. 7s = 1) [ZS]. While this model matches the STAR 
measurement at high energieqj, it does not provide a good description at the intermediate energies. On 
the other hand, if 7s is included in the fits as a free parameter, the yields can be described rather well 
(see right panel of Fig. [T7|) . This sensitivity of the strangeness neutral (j) meson to the 7s parameter 
suggests that the abundances of strange quarks are already determined in the partonic and not in the 
subsequent hadronic phase. 

In a hadronic scenario the can be produced via kaon coalescence. For instance in the hadronic 
transport model UrQMD |16] this is the dominant production mechanism. The predictions of this 
model are compared to the data as well (see Fig. [17]). At AGS energies the model under-predicts the 
(p/ir ratio at midrapidity by a factor of 2 (left panel) and is further below the measurements at SPS 
and RHIC energies. For the unnormalized total yields the discrepancy is smaller (right panel), but still 
prominent at higher energies. This leads to the conclusion that kaon coalescence is not sufficient to 
explain the amount of observed 0. There are two other observations that substantiate this point of 
view: the first is that the 0/K~ ratio does not increase with system size, as predicted by UrQMD (see 

"'Please note that there is an unresolved discrepancy between the STAR and PHENIX measurement. The latter is 
lower by 40% than STAR. The difference is essentially in the low pt region, while at higher transverse momenta the 
measurements are roughly in agreement. 
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lower left panel of Fig. [TK]) . but is essentially fiat. The second experimental indication is based on the 
comparison of the measured widths of the rapidity distributions to the widths expected in the case 
of kaon coalescence (see discussion in Section l4.3.2p . 

It is a remarkable observation that the 0/K~ ratio seems to be independent not only of the system 
size, but also of the center-of-mass energy of the reaction systems (see left panel of Fig. [TS]) |18U] . 
Alternatively, one finds that the midrapidity yields scale as well with the K"*" yield squared (middle 
right panel of Fig. [T8|l . Since the number of is directly proportional to the number of produced 
anti-strange quarks (= number of strange quarks), the (f) production thus might only depend on the 
product of s and s quarks. Since the s quarks are to a large extent contained in the K~ and A, the 
product s ■ s should be roughly proportional to K"*" ■ (K~ + A). The ratio 0/(K'*' x (K~ + A)) is similarly 
energy independent than 0/ (K"*")^ (lower right panel of Fig. [T51) and has approximately the same value 
(0.2 — 0.3 ■ 10"'^). On the other hand, the ratio 0/(K'''K~), as it might be suggested by a naive kaon 
coalescence scenario, is clearly not energy independent (upper right panel of Fig. ITSj) . 
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Figure 19: Left: The rrit spectra of Kg, A, A, S~, S"*", 0~, and fi"'" as measured by NA57 for central 
Pb+Pb collisions at 158 A GeV (Fig. adapted from |186] ). This solid lines represent the result of a fit 
with a hydrodynamically inspired model [187]. Right: The inverse slope parameters as a function of the 
particle mass measured in central Pb+Pb collisions at 40^4 GeV (closed symbols) and at 158A GeV(open 
symbols). The data are from the NA57 [HI] and the NA49 |59l [HH [IlO] collaboration. 



4.3 Spectra 

This section discusses the dependence of transverse mass spectra (mt = \/pt + ^^) on the center-of-mass 
energy of the reaction, which gives insight into the evolution of radial flow. This is of interest especially 
for multi-strange baryons, since they could reveal the onset of a partonic contribution. The inverse 
slope parameters extracted from mt spectra of charged kaons exhibit a step-like feature in their energy 
dependence, which might be connected to a first order phase transition. Further the measured rapidity 
spectra at different center-of-mass energies and attempts to connect thermal parameters determined at 
different rapidities and at different energies are discussed. 



4.3.1 Transverse momentum spectra 

Transverse mass spectra have been measured for all strange particle species at different center-of-mass 
energies. As an example, the left panel of Fig. fTU\ shows the mt spectra of Kg and hyperons, measured 
by the NA57 collaboration in Pb+Pb collisions at 158A GeV. Since strange particles span a relatively 
wide range of rest masses (between m(K^) = 493.677 MeV/c^ and m(f2) = 1672.45 MeV/c^), they 
provide an ideal probe for the effects of transverse expansion, which is generated by the high energy 
density in the fireball. Especially the question to which extent the radial flow is caused by the partonic 
stage of the reactions can be addressed with multi-strange particles. In order to quantify the spectral 
shape, the invariant mt spectra are usually fitted by an exponential function 

mt drritdy \ T* J 

with T* as the inverse slope parameter. For kaons the exponential results in a very good fit over a larger 
region of mt. Heavier particles exhibit deviation from an exponential behavior, due to the stronger 
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Figure 20: The transverse mass spectra of fi^ and fi^ in central Pb+Pb collisions at 158^4 GeV and of 
the sum fi^ + r2+ at 40^4 GeV [59j. The dashed curve shows a fit with an exponential. The solid, dotted, 
and dash-dotted curves represent a model including transverse expansion |187] . The used parameters 
are Tfo = 90 MeV and {l3±) = 0.5 (A), Tfo = 170 MeV and = 0.2 (B), and Tfo = 127 MeV and 
{I3±) = 0.5 (C). In (A) and (B) a linear velocity profile is used, while (C) was calculated with a constant 
expansion velocity. 

infiuence of radial fiow which introduces a convex shape in the spectra. Therefore, the fit result will 
depend on the rrit range in which the fit it performed. Nevertheless, a comparison of the inverse slope 
parameter T* extracted for different particles at different energies, reveals already interesting features. 
In the right panel of Fig. [12] such a comparison for 40^4 and 158^4 GeV beam energy is shown |188j . 
The inverse slope parameters increase linearly with particle mass up to the A. This can be understood 
as a consequence of the radial expansion of the fireball, which, in a simplified picture, introduces an 
additional component to the effective inverse slope parameter T*, that depends on the particle mass m 
and the expansion velocity v± as 



While the kinetic freeze-out temperature Tfo is assumed to be the same for all particles, the second 
term introduces the observed mass dependence of T*. The slightly higher T* measured at 158^4 GeV 
would thus indicate that the average radial expansion velocity still increases between the two energies. 
However, the S and Q do not fit into this systematics. At 158A GeV their inverse slope parameter is at 
the same level as the one measured for the A, while at 40^4 GeV it is significantly lower. Apparently, 
the multi-strange particles do not participate in the radial fiow to the same extent as the particles 
with one or no strange constituent quark. This also follows from the analysis shown in Fig. |2D] 
Here the rrit spectra of the Q are compared to a hydrodynamically inspired model, which assumes a 
transversely expanding emission source [187] . Such an analysis avoids the problems connected to the 
use of an inverse slope parameter, whose value for heavy particles can depend on the fit range in rrit. 
The parameters of this model are the kinetic freeze-out temperature Tfo and the transverse fiow velocity 
/3s = Vs/c at the surface. Assuming a linear radial velocity profile f3±{r) = f3s r/Rs, as motivated by 
hydrodynamical calculations, the rrit spectrum can be computed from 



T* ^ Tfo + -mvl 



(11) 




(12) 
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Figure 21: Left: The kinetic freeze-out temperature Tfo versus the transverse flow velocity for the 
hydrodynamically inspired model fits to the rrit spectra [132j. The 1, 2, and 3 sigma contours are shown. 
Solid curves are for a simultaneous fit to the S~ and Dashed curves are a separate fit to the vr, K, 
p, and A data as measured by the STAR collaboration. The diamonds represent the best fit in both 
cases (Fig. adapted from |132j ). Right: The mean transverse momentum {p±) for identified particles 
as a function of the particle mass |132j . The band follows from the three sigma contour of the fit to tt, 
K, p, and A, as shown in the left panel, and the dashed curve is for Tfo = 170 MeV and {(3±) = 0. 

where Rs is the radius of the source and p = tanh.~^(3± the boost angle. The curve labeled (A) 
is calculated with the parameters derived from a simultaneous fit to kaons, (anti-)protons A, and 
(Tfo = 90 MeV and {f3±) = 0.5). The clear disagreement to the measured Q rrit spectra illustrates that 
the freeze-out conditions for the VL are different than for the lighter hadrons. A much better agreement 
can be achieved, if the freeze-out parameters taken from a fit to ^/ip and ip' spectra |191j are used, 
as shown as curve (B) in Fig. [201 This indicates that the VL has similar kinetic freeze-out conditions 
than the J/V"- The same observation was made at RHIC energies (see Fig. Using the same model 
different freeze-out parameters are extracted for the H than for the lighter hadrons. The averaged 
transverse momentum of the S and VL departs from the linear increase with particle mass, similar to 
the inverse slope parameters shown in the right panel of Fig. [191 

The interpretation of this behavior follows from the assumption that rare particles as the S and 
have a lower hadronic scattering cross section than light hadrons |192j and therefore do not participate 
in the radial flow that is developing during the hadronic phase of the flreball evolution. This leads to 
the conclusion that a substantial part of the transverse expansion probed by these particles has to be 
generated during the partonic phase. Thus, the S and VL would be directly sensitive to the pressure in 
the early phase of the reaction. The rapid increase of the inverse slope parameters measured for these 
particles in the SPS energy range (right panel of Fig. [T9|) would thus indicate that the partonic part of 
the flow is beginning to develop at these energies. 

The mt spectra of kaons provide important information on the transverse dynamics of the reaction 
system. Calculations within hybrid hydrodynamics and transport approaches [1931 1194] indicate that 
kaons freeze out slightly earlier than nucleons and A (e.g. in |195] decoupling times of Tfo ~ 14 fm/c are 
found for mesons and Tfo ~ 18 fm/c for nucleons, A and E), and thus they will be less sensitive to the 
late rescattering phase. Their mt spectra can be very well described by a single exponential (Eq. ([TO|l ) 
and can thus be characterized by just one parameter, the inverse slope parameter T*. The energy 
dependence of T* for charged kaons exhibits an interesting feature, as shown in Fig. [221 While T* is 
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Figure 22: The energy dependence of the inverse slope parameters T* of K"*" (left) and K~ (right) 
measured at midrapidity in central Pb+Pb and Au+Au collisions |126] . The K"*" slope parameters are 
compared to results for p + p(p) collisions |196] (left). The curves shown in the right panel represent 
various model predictions. 

rising rapidly with center-of-mass energy for y/s^ < 7 — 8 GeV, it is rather constant or only slightly 
increasing above this energy. A similar observation has been made for the averaged transverse mass 
(mt) — m of pions and protons |126[ 1197] . In proton-proton collisions such a behavior is not observed 
(see left panel of Fig. [22]) [T96] . 

The dN/dy dependence of the mean transverse momentum has been suggested as a signature for 
a first order phase transition to a deconfined state quite a while ago within the concept of Landau's 
hydrodynamical model |198] . Originally, this has been proposed as an interpretation of the dependence 
of (pt) on dN/dy observed in p + p collisions at y/s^ = 540 GeV. While the fireball volume produced 
in these reactions is most likely too small to create a QGP state, attempts have been made to apply a 
similar interpretation to the observed structure in the energy dependence of the inverse slope param- 
eters of kaons [199J. While hadronic transport model calculations are generally not able to properly 
reproduce the observed structure (see right panel of Fig. |22]) |l6l [55 | 1168] . the data can be described by 
a hydrodynamical model assuming a first order phase transition between a hadronic and a deconfined 
phase (dash-dotted curve in right panel of Fig. |22|) |199j . If initial conditions calculated event-by-event 
with the NEXUS event generator [2001 1201] are used, the same model is able to simultaneously describe 
yields, transverse mass and rapidity spectra. Another hydro approach, which uses a dynamical descrip- 
tion of the freeze-out conditions [202] . also results in a good match to the data, but turns out to be 
less sensitive to a specific equation-of-state and thus does not necessarily require a phase transition. A 
similar observation has been made with a model that combines hadronic transport and hydrodynamics 
|203] . Also here it was found that the way the freeze-out procedure is implemented can have as much 
influence on the final result as the equation-of-state that is used. Even though it is quite possible that 
the measured energy dependence of the inverse slope parameters might be due to a first order phase 
transition, it is currently difficult to establish an unambiguous theoretical connection. At the moment 
it is also still unclear whether in the center-of-mass energy range discussed here a first order phase 
transition is to be expected. While most lattice QCD calculation expect a change from a cross over to 
a first order phase boundary at higher /Ib, corresponding to lower ^ys^^ [21 IH IS], the results obtained 
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Figure 23: The rapidity spectra of vr", K^, K~, 0, A, and A as measured by the NA49 collaboration in 
central Pb+Pb collisions at 158 A GeV (197j . The closed circles indicate measured data points, while 
open ones are reflected around midrapidity. The solid lines represent fits with a single Gaussian or the 
sum of two Gaussians. 



by [B] rather indicate that the phase transition is of the type of a cross over also at higher /Xb- 



4.3.2 Rapidity spectra 

At the SPS rapidity spectra in the forward hemisphere have been measured by the NA49 collaboration 
for charged kaons, 0, and hyperons at five different beam energies [581 EHl 11251 I126i 1149] . Part of 
these spectra are shown for central Pb+Pb collisions in Fig. [231 The NA57 collaboration performed 
measurements of rapidity distributions for Kg, A, A, S^, H^, and ^2"+^"*' in the lab rapidity range 2.5 < 
y < 3.5 for the 53% most central Pb+Pb collisions at 158yl GeV |2U4j . At RHIC no measured rapidity 
distributions are available for most particle types, since most detectors only cover the midrapidity 
region. The only exception is provided by the BRAHMS collaboration, which can cover a large part of 
the forward hemisphere. BRAHMS has measured in central Au+Au collisions at y/s^ = 64.2 GeV and 
200 GeV rapidity distributions of dN/dy and (pt) for charged kaons (see Fig. |2l]) |166[ 1174] . In order to 
avoid assumptions on longitudinal phase space distributions, the measurement of dN/dy distributions 
over a large range of rapidities is a prerequisite for the precise determination of Att integrated total 
yields. Therefore, measurements of the total multiplicities at RHIC energies are only available by the 
BRAHMS collaboration for charged pions and kaons (see Tables [3]- [TT|) . 

The rapidity spectra of strange mesons can reasonably be described by a Gaussian distribution at 
all energies |149[ 1166] (Figs. |23l and [26]) , although in some cases the sum of two Gaussians, displaced 
relative to midrapidity, can result in a slightly better fit |125[ 1126] . This observation is also true for 
anti-baryons (A, S"*") [SB] and for the Q [59j. Baryons, such as the A and S~, have a wider distribution 
which is changing drastically with beam energy, because their light valence quark structure is connected 
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Figure 24: dN/dy (upper panels) and the averaged transverse momentum (pt) (lower panels) as a 
function of rapidity for charged pions and kaons for central Au+Au collisions, as measured by the 
BRAHMS collaboration at = 62.4 GeV (left) [HI] and at = 200 GeV (right) [T66] . 

to the ones of the originally colliding nucleons. This is in particular visible for the A, whose distribution 
changes from a Gaussian shape at 20A GeV to a spectrum that is essentially flat inside the measured 
region at 158 A GeV (see lowest row of Fig. 1251) . No pronounced change with center-of-mass energy is 
observed for S~ spectral shapes, but the widths of the rapidity distributions increase significantly faster 
with y/s^ than the ones measured for S"*" [58j. The reason for this behavior is that A and S~ are 
sensitive to the longitudinal redistribution of net-baryon number, which depends strongly on center-of- 
mass energies. While at low SPS energies the net-baryon number is concentrated around midrapidity, 
its distribution changes very rapidly and turns from a maximum around y = into a shallow minimum 
at RHIC energies |205[ 1206] . Consequently, the midrapidity anti-baryon/baryon ratios exhibit a strong 
energy dependence. With increasing strangeness content of the baryons this effect becomes weaker, as 
already shown in the right panel of Fig. |9l 

The systematics of the RMS widths of the rapidity spectra of vr^, K^, 0, and A, i.e. those particles 
whose longitudinal distributions can be described by Gaussians, is summarized in the left panel of 
Fig. [25] |149] . For charged pions, kaons, and anti- lambdas the widths a scale linearly with the beam 
rapidity ybcam = 2 cosh~^(ys^/(2 mp)) in the SPS energy region |197] . The widths depend, at least 
in the SPS energy region, on the particle type as cr(7r~) > cr(K"^) > cr(K~) > cr(A). However, since 
at RHIC the width of the K"^ approaches the one of the pions, or is even larger jl66] . this scaling 
apparently does not hold up to top RHIC energy. 

A peculiar behavior can be observed for the meson, whose rapidity widths are increasing much 
faster with beam energy than the ones of the other mesons (left panel of Fig. [2^ . While at 20^4 GeV 
the has a rapidity distribution that is narrower than the one of the K~, its width is comparable to 
the one of the vr"*" at 158A GeV. This behavior is at variance with a simple expectation based on kaon 
coalescence (right panel of Fig. [25]). In this case the rapidity distributions a^f, should depend on the 
measured widths of the kaon rapidity distributions aK± in the following way: 

111 
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Figure 25: Left: The RMS widths of the rapidity distributions of 7r~, K~, K , and (j) in central Pb+Pb 
colhsions at SPS energies as a function of the beam rapidity |125[ I126[ 11491 1197] . Right: The widths of 
the rapidity distributions of mesons in central Pb+Pb collisions at SPS energies as a function of the 
beam rapidity |149] . The data are compared to expectations from kaon coalescence (solid line). (The 
Figs, are adapted from |149] .) 



As shown in the right panel of Fig. |25l this simple scenario is clearly ruled out for nucleus-nucleus 
collisions. However, it is interesting to note that for p+p collisions at top SPS energy (158 GeV) 
the coalescence expectation, derived from the measured kaon rapidity widths according to Eq. f lT3|) . is 
relatively close to the measured rapidity widths (a™^'^^' = 0.89 ± 0.06 |207] . compared to cr^"^^' = 0.77 
as based on the kaon measurements in |208j ). In addition, it is remarkable that the rapidity widths 
measured at top SPS energy are significantly wider in central Pb+Pb collisions compared to p+p 
collisions at the same energy (/2MS'y (Pb+Pb) = 1.44+0.021+0.054 [HH] and i?MS'y(p+p) = 0.89+0.06 
|207] . see also left panel of Fig. [25]). For instance, such a significant difference between p+p and A+A is 
not seen for charged kaons (K+: /2MS'y(p+p) = 1.20 and i?MS'y(Pb+Pb) = 1.22, K": i?MS'y(p+p) = 
1.01 and i?MS'y(Pb+Pb) = 1.14, both at = 17.3 GeV [1251 M)- 

In a publication by the BRAHMS collaboration |166] the K/tt ratios measured at different rapidities 
for central Au+Au collisions at y/s^ = 200 GeV were compared to the corresponding ratios determined 
at midrapidity, but at different center-of-mass energies. While at top RHIC energy the K~ /tt~ ratio 
drops when going from midrapidity to forward rapidities, the K'^ /7r~^ ratio rises. Since the production 
is favored by the presence of baryons via associated production mechanisms (e.g. p + p — )■ p + K"*" + A), 
the higher net-baryon number at forward rapidities at RHIC can thus increase the number of relative 
to K~ and explain the rise in the K^/tt"'' ratio. As the net-baryon number at midrapidity is increasing 
when going toward lower center-of-mass energies, the K^^/tt^ ratio should rise to a similar level when 
the same net-baryon number is reached. Figure 126) contrasts the two dependences. The K~^/7r^ ratio 
measured at very forward rapidities for ^,/s^ = 200 GeV roughly corresponds to the one observed at 
midrapidity and y/s^ = 8.7 GeV. There seems to be a universal dependence of the K~/K^ratio on the 
p/p ratio, if measurements at different energies and rapidities are compared (see Fig. [T6|l . 

In the context of statistical models, this dependence on the net-baryon density is described by 
the baryo-chemical potential /ig. If the chemical freeze-out temperature is assumed to be constant, 
which is a good approximation for y/s^ > 17.3 GeV , the change in /ig caused by the changing net- 
baryon number or the changing center-of-mass energy would therefore be the only parameter driving 
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Figure 26: The K/vr ratios, measured in full phase space for central Pb+Pb/Au+Au collisions, as a 
function of ^ys^ (left) and of rapidity (right) at y/s^ = 200 GeV |166] . The dotted and dashed lines 
correspond to predictions from a statistical model with rapidity independent freeze-out parameters 
[209] . 

the changes in the other particles ratios. However, this is not a universal behavior valid everywhere. 
For instance, the /n'^ ratio is essentially independent of rapidity at y/s^ = 62.4 GeV, even though 
the net-baryons are changing quite dramatically with rapidity |174] . Assuming that there are no large 
differences between the vr"*" and vr" rapidity spectra (only vr" spectra have been measured), there are 
indications that for Pb+Pb collisions at ^/s^ = 17.3 GeV the /Tr~^ ratio drops instead of rises with 
increasing rapidity, which is opposite to the behavior seen at y/s^ = 200 GeV. This can be deduced 
from the fact that at y/s^ = 17.3 GeV the K"*" rapidity width is smaller that the one for pions (see 
left panel of Fig. [2^ ^125] , while it is wider at top RHIC energy (see right panel of Fig. [21]). Since 
at y/s^ = 17.3 GeV the net-baryon number is slightly rising toward forward rapidities (at least up 
to = 1.8 |205j ). this would indicate that it is not the only driving factor for the relative strangeness 
production. As already mentioned in Section I4.2.2[ a study of the rapidity dependence of the chemical 
freeze-out parameters derived from fits with a statistical model shows that both, Tch and /^b, do depend 
on rapidity |179] . In the same analysis, it has been found that the dependences on y are not the same 
at different center-of-mass energies. The observed dependence of statistical model fit results on whether 
midrapidity dN/ dy or 47r integrated data are used as input (see discussion in Section 12. 3p seems to be 
due to this effect. Strangeness production in the fragmentation region is apparently different than in 
the central fireball and it starts to affect the measured rapidity distributions for y/s^ below top SPS 
energies. 
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5 System size dependence 



Strange particle production rates in heavy ion collisions, relative to elementary p+p collisions define 
the strangeness enhancement (Eq. (Q) as already discussed in Section \3A\ of this article. Whether 
the strangeness enhancement in larger systems is due only to the phase space suppression in the smaller 
systems, such as p+p or p+A, or is indeed due to an increased strangeness production from a deconfined 
partonic medium is still an ongoing discussion. 

An rapid increase of the strangeness enhancement factor in small systems (up to A^part = 50) has been 
suggested by canonical models due to the restriction in phase space. In other words the volume is not 
sufficiently large to contain the necessary strangeness content for statistical strange particle production. 
Complete statistical production can be assumed when the system reaches its chemical equilibrium and 
is thermalized. Strangeness equilibration requires a finite number of binary collisions in an extended 
volume. This can only be achieved if the number of participants is larger than 50. In comparison, 
the light quark hadrons equilibrate faster, i.e. in a smaller system with less binary collisions. In p+p 
collisions the underproduction of strangeness due the limited number of binary collisions (one binary 
collision) is largest, which means that strangeness is suppressed, at least for collision energies up to 
RHIC energies. 

Any centrality scaling of strangeness production with respect to the p+p yield therefore includes a 
strangeness suppression factor. This means that any apparent strangeness enhancement factor Es as 
a function of centrality is foremost dominated by the transition of a non-equilibrated small heavy ion 
system into an equilibrated large system. Any remaining additional enhancement could be interpreted 
as a phase transition signature. The strangeness enhancement as a function of system size, A^part and 
as a function of collision energy is shown in Fig. [27] (left) for the fl in Pb+Pb collisions |210] . 
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Figure 27: Left: The enhancement of inclusive protons. A, A, S~, S"*", and Q~ + 0,^ around midrapidity 
as measured by the STAR collaboration for Au+Au colhsions at ^/s^ = 200 GeV (solid symbols) |119] 
and by the NA57 collaboration for Pb+Pb collisions at ^/s^ =17.3 GeV (empty symbols) |112] . The 
boxes at unity show the combined statistical and systematical errors in the p+p, respectively p+Be, 
reference data. The arrows mark statistical model predictions for a grand-canonical ensemble and two 
different temperatures (T^h = 165 MeV and Tch = 170 MeV). Red arrows correspond to S~ (S+), 
black arrows to A (A). Right: Centrality dependence of relative enhancement of the yields of Q~ + Cl^ 
normalized by the number of participants in central Pb+Pb to p+p reactions at different collision 
energies |210] . 
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When the enhancement factor levels out it suggests that the collision region has reached sufficient size 
and the small volume size effects can be neglected. In statistical models, the predicted equilibration of 
the strangeness production, independent of the incident energy, is realized for sufficiently large volumes 
(-^part > 50). At the highest SPS energy and the RHIC energies the Q yield increases by about an order 
of magnitude due to this phase space effect alone. Lower energies yield an even larger enhancement 
factor due to the more dominant suppression of strangeness production in small systems. Any additional 
enhancement due to deconfinement would have to be assessed on top of this 'trivial effect', i.e. it should 
manifest itself as an enhancement relative to the equilibrium value. Since it is assumed that at RHIC 
energies and below the QGP will arise only in sufficiently large systems, a systematic study of the volume 
dependence of the enhancement factor might give information about the onset of QGP creation. 

Recent results from the LHC in elementary collisions might indicate that in high multiplicity events 
a phase transition might occur in p+p collisions. At RHIC energy, though, this evidence could not be 
obtained due to the low collisions energy |211l 12121 1213j . 

The dependence of strangeness production on the system size (volume) of the collisions has been 
investigated via two methods, either by the variation of the target and beam nucleus species or by the 
selection of the centrality (impact parameter) of the collisions. In order to characterize the collision 
geometrical calculations such as Glauber |115[ 1116] models are used to determine the interaction volume, 
the number of participating nucleons ( = wounded nucleons) and the number of binary collisions in the 
heavy ion reaction. 

Measurements from STAR and NA57 at y/s^ = 200 GeV and y/s^ = 17.3 GeV are shown in left 
panel of Fig. 1571 They reveal a steady increase of strange hadron production with increasing system 
size (iVpart) and increasing strange quark content. A saturation (leveling off of the enhancement factor) 
for large systems as predicted by statistical models is not detected. In fact, the enhancement factor 
continues to increase up to the most central collision (A/'part = 350). This observation can have different 
explanations. For instance, it might suggest that even in the most central events the full equilibrium 
value of a grand-canonical ensemble is not yet reached for multi-strange particles. 
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Figure 28: Temperature (left) and strangeness saturation factor (7s) (right) as a function of the number 
of participants (A^part) for all particles measured by the STAR experiment |214l 1216] . as obtained with 
the THERMUS code [63. 
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An alternative explanation for the enhancement above A^part = 50 would be the observation of 
strangeness enhancement due to a QGP phase in a grand-canonical system, i.e. in systems where the 
volume effects of the strangeness suppression from p+p is not present anymore. Calculations of the 
strangeness saturation factor 7s (Fig. 128) right), using a thermal model fit to STAR data |214[ 1215] are 
consistent with this interpretation. The saturation factor 7s approaches unity for Apart = 100, which can 
be interpreted as evidence for an equilibrated and strangeness saturated system at these system sizes. 
Therefore, any further enhancement of strangeness production from Apart = 100 on might be related to 
an additional contribution from a deconfined medium. At RHIC energies this could account for about 
half of the enhancement factor for Q baryons measured in central Au+Au collisions and compared to 
p+p collisions. 

5.1 Core and corona model 

A different conclusion can be reached if one considers a description of the interaction region in a heavy 
ion collision, which distinguishes between independent core and corona volumes in centrality dependent 
heavy ion collisions |157[ I217[ I218[ 1219] . Here the dense core behaves like a QGP volume, whereas 
the surrounding corona can be understood as a superposition of independent p+p collisions. The 
justification for such an approach is based on the realistic modeling of the nuclear surface rather than 
the very schematic rigid sphere calculations used in earlier models. 

The effect of core-corona models is best investigated in the centrality dependence of particle pro- 
duction since the corona fraction continues to increase for the more peripheral systems. In other words, 
the core fraction, which represents the QGP volume, increases with centrality. The system might be 
equilibrated in a grand-canonical regime even at small Apart but due to the volume to surface ratio 
increase the strangeness enhancement factor continues to increase. 

In order to investigate the strangeness enhancement, this volume effect of the core-corona ratio 
has to be taken into account and subtracted from the measured strangeness enhancement factor. The 
remaining enhancement will be the real "strangeness enhancement". 

5.1.1 Centrality dependence 

To investigate the strangeness enhancement as a function of the collision centrality the additional effects 
of the change in surface to volume ratio will be taken into account by unfolding the source into a corona 
(surface) and core term. The total volume is the superposition of the core and corona volume. Both 
areas will be treated independently, i.e. there are no interactions between core and corona particles. 

As a consequence of this superposition in heavy ion collisions one has to calculate the relative weight, 
in terms of particle production, of core and corona separately as a function of centrality. The particle 
production of the core volume is described by a partonic medium in thermal equilibrium whereas the 
corona is described by the nucleon-nucleon collisions modeled by a statistical hadronization model |157] 
(Becattini) or alternatively by the EPOS model |218] . 

Figure 129) (left) shows the corona participants (Apc) as a function of the number of total participants 
Apart, which increases sightly with centrality. Although the number of corona participants increases their 
relative fraction of (Apc/-^part) drops dramatically from around 50% at A'part = 50 to 10% at Apart = 350. 
This fraction can be calculated either with a Glauber model (red squares) or from a fit (green points) to 
the strange particle yields of K, A, 0, S and Q (Fig. [3T1) when using relative contributions from corona 
and core type mechanisms as shown in Fig. [29] (right). The fitted A^pc and the A^pc from Glauber 
model calculations are in agreement. However the errors on the corona participants A^pc derived from 
the data are very large and therefore not sensitive to falsify this approach. 

To quantify the relative contribution of the two components, the fraction of nucleons that scatter 
more than once /(A^part) can be used. /(A^part) can simply be calculated with a Glauber model |115] . 
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Figure 29: Left: Number of corona participants A^pc at different centralities as a function of the total 
number of participants A^p. The square dots denote the values calculated with a Glauber model while 
the round ones denote the values arising from fitting A^pc as a free parameter. The round symbols are 
shifted 5 units of A^p leftward for clarity. Right: Fraction of produced particles coming from core (upper 
lines) and corona (lower lines) as a function of centrality |220] (Becattini). 




Figure 30: Multiplicity per participant as a function of A'^p for the core part only. The results for tt , 
K~, p, A, S+ in Au+Au collisions at ^/s^ = 200 GeV (dotted lines), and for n~, K~, Kg, A, S+ in 
Pb+Pb collisions at = 17.3 GeV (dashed lines) (Fig. is adapted from |2lH] (EPOS)). 



This quantity allows for a natural interpolation between the yields Y measured in elementary p+p 
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(= ^corona) and in central nucleus- nucleus collisions (= Icore)- 

Y{Np,,t) = iVpart [/(iVpart) l^core + (1 - /(iVpart)) Vcorona] (14) 

The contribution of the core part alone is shown in Fig. [30] for a calculation using the EPOS model, 
which can be compared to the result from an alternative model by Becattini shown in Fig. 121] (right). The 
overall trend of increased particle production from the core volume of this models is the same. However 
the particle yields differ, mostly in the region of small number of participants. This discrepancy can be 
investigated experimentally after statistically more significant data samples are obtained in particular 
for either small collisions systems or peripheral collisions of the heavy systems. 

In summary, the derived strangeness enhancement factors from different core-corona models shown 
in Fig. EI] (left: Becattini, right: EPOS) describe the steady increase of the data from p+p up to the 
most central heavy ion collisions. 

The S+ yields in the more central collisions are better described by Becattini than the EPOS Model. 
However, this model is over-predicting the S"*" yield for peripheral collisions in the range of A^part < 60, 
which might be a hint that the centrality dependence cannot exclusively be described by a core/corona 
volume effect. 

A final interpretation of this result will depend strongly on the modeling of the core volume. If the 
dynamics of the produced particles require the system to be partonic in nature, which is also indicated 
by measurements of the radial and anisotropic expansion of the system, then the enhanced strangeness 
production is not just a volume effect, but indeed a signature for deconfinement, at least in the core 
of the collision volume. The main contribution of the core-corona models to this discussion is a more 
precise understanding of the centrality dependence of the strange particle production which goes beyond 
the simple phase space argument of strangeness suppression in small systems. 
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Figure 31: Left: Hyperon rapidity density measured in Au+Au collisions normalized to measurements 
in p+p collisions and the number of participants at y/s^ = 200 GeV as a function of the number 
of participants |119j . The lines are theoretical calculations from Becattini et al. |157[ 1220] . Right: 
Rapidity density dN/dy per participant as a function of the number of participants (iVp) in Au+Au 
collisions at y/s^ = 200 GeV (RHIC) for tt~ , K~, p. A, S+ together with the full calculation (full lines) 
and the core part only (dotted lines) (Fig. is adapted from |218] (EPOS)). 
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5.1.2 Comparison of different collision systems 

A measurement which might also be explained by the core-corona approach is the apparent difference 
of the strangeness enhancement factors for different collision systems at the same iVpart- This was first 
measured at the SPS |221j and has since then been corroborated at RHIC energies [222]. 
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Figure 32: The enhancement factor for (multi-)strange particles at midrapidity in Cu+Cu and Au+Au 
collisions at y/s^ = 200 GeV. The A and A yields have been feed-down subtracted in all cases. The 
Au+Au multi-strange data have been previously published jl33] . The black bars show the normalization 
uncertainties, and the uncertainties for the heavy ion points are the combined statistical and systematic 
[222] . The lines represents calculations with the EPOS model |218] . 
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Figure [32] shows the difference between the Cu+Cu and Au+Au systems at y/s^ = 200 GeV colhsion 
energy. For all strange particle species the smaller system exhibits larger strangeness enhancement at 
the same A^part SiS is summarized in Fig. 1321 A possible explanation is the larger core in Cu+Cu compared 
to Au+Au when using a core-corona approach. At the same Apart the Cu+Cu collision system features 
a larger core fraction and thus a larger number of participants that undergo multiple collisions according 
to Glauber calculations. This is implemented in EPOS which describes the trend of the data in Cu+Cu 
better than in Au+Au collision in Fig. [32j On the other hand the effect is also visible in the AMPT 
model which does not feature separate core and corona volumes but rather uses string breaking and 
mini-jet fragmentation to describe the medium. In general though it seems true that at the same 
A^part the Cu+Cu system yields a larger contribution to the strangeness yield than the more peripheral 
Au+Au system. 



5.2 Scaling with A^part and A^bin 

In the case of a chemically equilibrated bulk system of light quarks we expect rather A^part scaling 
than A^bin scaling for the centrality dependence of the particle yields. This is true, for example, for 
the protons at RHIC energies (Fig. [27] left). Their enhancement factor, which presumes A^part scaling, 
is constant as a function of the system size. In other words, the increase in proton production from 
peripheral to central heavy ion collisions scales with the number of participants. This hints at the 
dominant production of up and down quarks through thermal parton-parton interactions rather than 
hard scattering. 

A^part scaling does not necessarily proof that the system in its early stage is thermalized. A particle 
species, whose initially produced yield scales with A^bm and is thus not chemically equilibrated, might 
thermally equilibrate in the hadron phase through elastic scattering and exhibit flow features and a 
thermal pt spectrum in central collisions. On the other hand, species for which the yield scales with 
A^part might exhibit non-thermal behavior in their kinematic spectra. 

Generally heavier quarks (i.e. charm and bottom) are not expected to equilibrate chemically or 
thermally, and are mostly produced through initial parton-parton interactions which scale with the 
number of binary collisions, A^bin- Since the strange particle yield continues to increase beyond the 
Apart scaling relative to p+p, this could be viewed as a non-equilibrium signature for strange quarks. 
An alternative scaling method, based on an independent normalization of the constituent quarks in a 
hadron was suggested by H. Gaines. The method takes into account the specific scaling behavior of the 
strange quark, which is the A^bin scaling. The new factor is defined as |223] : 

^ _ A^light X A"part , A"s X A"bin , . 

Oscalmg- + libj 

where Aq is the number of quarks in the particle, A^ught is the number of light (u and d) quarks, and 
As is the number of strange quarks. Figure [33] right shows the relation of Apart and Abin for different 
collision energies based on a Glauber calculation [1151 1116] . This numbers are used to calculate the 
new enhancement factors in Fig. [33] left, which shows a fiat distribution for all light and strange quark 
particle species. 

This new scahng is successful on the level of 20%. Thus one could conclude that the strange quark 
indeed scales with Abin and is thus mostly produced through initial gluon-gluon interactions. However, 
one exception is the scaling of the 0, which according to Eq. ([TBI) should scale with A'bin, but appears 
to scale better with Apart- 
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Figure 33: Left: Au+Au data scaled as defined by Eq. f|T5l) except the 0, which is scaled by A'part- 
All data are normalized to the most central bin [223] . Right: Number of participants versus number 
of binary collisions using Glauber calculation for Pb+Pb and Au+Au collision systems at Js^^ = 
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200 GeV [miTTB] . 



6 Strangeness at intermediate and high 

The bulk of the particles (> 90%) produced in heavy ion collisions carry low momentum {pt < 2 GeV/c) 
and are produced either through soft multiple parton scattering or statistical production from a thermal 
medium. The intermediate pt range is commonly defined as the region between pt = 2 and 6 GeV/c 
and can be considered the transition region between hard processes, described by pQCD, and the 
aforementioned soft region. Hadron production from pure fragmentation of the hard scattered parton is 
only dominant at pt > 6 GeV/c based on the RHIC heavy ion results. In a medium we expect a modified 
parton fragmentation 1225^226] 1227] which presumably softens the fragmentation function, i.e more soft 
particles are produced in the emitted hadron jet and pushed into a region which is dominated by the 
physics of the bulk. In this picture the particles in the intermediate pt range can be due to either a large 
radial expansion which adds additional kinetic energy to the bulk or an enhanced quenching of high 
Pt jet fragments traversing the medium. An additional alternative arises if one takes into account the 
possibility of forming color neutral objects through recombination of partons in the deconfined medium. 
This recombination mechanism for hadronization has been first introduced for elementary collisions in 
the late 70's [28j. When applied to heavy ion collisions the initial system is a thermalized partonic 
medium in which partons close in phase space can recombine to form a hadron with momentum from 
the sum of the quark momentum. 

Since the thermal spectrum of partons is likely to have a slightly higher mean momentum than the 
pions from the kinetic freeze-out, and the total momentum of the final hadrons depends on the number 
of added quarks, this mechanism populates predominantly the mid-pt range and favors baryon over 
meson production at a fixed intermediate pt [29l [30]. Thus the two main experimental measures to 
quantify the interplay of bulk, fragmentation and recombination processes are the pt dependent baryon 
to meson ratio and the pt dependent nuclear suppression factors. Here we will first discuss the baryon 
to meson ratios. 
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6.1 Baryon to meson ratios 



The enhanced baryon/meson ratio in the transition from pt = 2 — 6 GeV/c in heavy ion colhsions 
indicates that particle production from the deconfined medium is neither simply statistical in nature 
nor can it be described through fragmentation alone. The earliest RHIC measurements of the p/vr ratio 
show an increase up to pt = 3 GeV/c to unity |127l 1228] . a value much higher than in elementary p+p 
collisions [2291I230] . 
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Figure 34: The A/Kg ratio as a function of pt for different collision centralities for Au+Au collisions at 
^ = 200 GeV [221]. 

For the strange sector Fig. shows the A/Kg ratio in Au+Au collisions at ^/s^ = 200 GeV out 
to Pt = 6 GeV/c for different centralities |224] . The maximum of the distribution is at pt = 2.8 GeV/c 
for all collision systems including the elementary p+p collisions. However, the maximum value of the 
ratio for p+p collision is below unity, while it continues to increase for more central Au+Au collisions 
up to 1.8 in the 5% most central collision. This remarkable value means that in certain momentum 
regions baryon production dominates over meson production by almost a factor of two. Such an effect is 
impossible to model in fragmentation where the baryon production relies on the formation of di-quark 
out of the fragmenting proton and is thus considerably reduced in probability. 

Early attempts to describe this phenomenon tried to merge a radially expanding bulk system, which 
pushes higher mass particles (baryons) to larger momenta than lower mass particles (mesons), with 
a common quenching model for high momentum hadrons |231] . Such a model qualitatively leads to 
a 'pile-up' of the baryon to meson ratios at mid-pt, although a quantitative comparison has never 
been attempted. On the other hand thermal recombination models became more quantitative over the 
years and now also take into account unique strangeness aspects in these particle ratios. The general 
recombination approach is based on the idea that since the thermal parton spectrum is exponentially 
falling as a function of momentum, the probability to form a high momentum three quark state is 
higher than the forming a two quark state of the same final state momentum. The effect is per se not 
unique to strange particles as shown by p/vr ratio, but certain aspects of the process can be tested 
when comparing to non-strange ratios. First, the models require thermalization of the strange parton 
in the system prior to hadron formation [31]. Second, the relevant degree of freedom in those models 
is the constituent quark, i.e. the modeling requires quarks to have a finite mass which differs between 
u, d and s-quarks by about 100 MeV/c^ [32]. More recent models even include recombination between 
bulk partons and the non-thermal hard partons ^B] . These models are successful in describing the pure 
strange ratio of f2/0 up to pt ~ 4 GeV/c as shown in Fig. |35] [232] . which includes more than 95% 
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of the hadron yield. This is of particular relevance since the small hadronic production cross section 
for Q and particles requires the particle generation to occur through partonic rather than hadronic 
recombination. At higher pt the model fails for the coalescence of thermal quarks and the combination 
of thermal and shower (total) quarks. Only the decreasing trend of from the total quark (thermal 
and shower) combination is visible. Since most of the 0-mesons are made via coalescence of seemingly 
thermalized s-quarks in central Au+Au collisions, the observations imply that a hot and dense matter 
with partonic collectivity has been formed at RHIC energies. 

Although coalescence models assume particle production from a source of thermal partons their 
results can be incompatible with thermal equilibrium of the final hadrons in the low (bulk) region. 
They should therefore only be applied to hadrons in a medium pt range (e.g. 2-6 GeV/c). In order 
to not violate entropy conservation, another production mechanism is required to produce a thermal 
reservoir of bulk particles below 2 GeV/c. Above ~ 6 GeV/c it is likely that parton fragmentation will 
be the most dominant mechanism. 
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Figure 35: The N(fi)/N(0) ratio vs. pt for three centrality bins in y/s^ = 200 GeV Au+Au collisions. 
The solid and dashed lines represent recombination model predictions for central collisions [33] for total 
and thermal contributions, respectively [232]. However the coalescence model should be applied for 
Pt> 2 GeV/c only. 



6.2 Nuclear suppression factors {Raa and Rcp) 

In order to study the medium modification of the hard process one has to determine the yield as a 
function of momentum. These pt dependent studies can be quantified via the nuclear modification 
factor defined by the yield in p+p {d^ a^^ / dptdy) and A+A (d^A^^^/dptdy) collisions, normalized by 
the binary collision scaling factor A^bin and the inelastic cross section for p+p collisions (0";^^): 

_ d'N^^/dp,dy 
"^^^ - ATAA dVNN/dp,d|/ ^'^^ 

Since the accumulation of sufficient statistics in the elementary p+p collisions required several run 
periods, the early RHIC results on nuclear suppression were based on Rcp, which is the momentum 
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dependent yield of central {d'^ N^^^^ / dp^dy) divided by the peripheral (d^A^P^^'/dptdy) collisions and 
normalized by their respective number of binary collisions (Nun) (see Fig. [361 right). 



_ NZ d^iV-Vdptdi/ 

AT-f d^ATP^/dptdy ^ ' 

For non-strange light quark particles this observable is roughly equivalent to -Raa since the light 
flavor production scales well from p+p to peripheral A+A collisions. In the strange sector, though, one 
would expect a distinct difference between Rcp and Raa based on the strangeness suppression in the 
elementary collisions which was discussed in the previous chapter. Rqp attempts to eliminate this phase 
space effect by normalizing the central spectrum to a more peripheral spectrum where one expects the 
strangeness suppression to be already minimized. 




Figure 36: Left: Raa for p+p |233l|23l], A and for a) 0-5% and b) 60-80% most central Au+Au 

events at ^/s^^ = 200 GeV. Right: Rcp measured in Au+Au events for p (0-12%/60-80%) [231], A+A 
and S-+H+ (0-5%/60-80%) [133] • The dashed curve shows the results for h++h- for 0-5%/60-80% 
|235] . In all plots the statistical and systematical errors are included. The band at unity shows the 
uncertainty on the binary scaling A'bin and the dashed line indicates the participant scaling (iVpart) with 
its uncertainty in a gray band. 
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Fig. |3S]right shows that the -Rcp of ^ is in agreement with the Rqp of protons out to pt = 3.5 GeV/c, 
while the Rqp of the S shows shght enhancement in the higher pt region. However, the errors are too 
large to draw any conclusion about an enhancement of strange particles with respect to non-strange 
particles. 

Fig. I3S] left shows -Raa for A and S compared to inclusive p+p measurements [2331 1234] . The 
Raa distributions of strange particles reach a maximum greater than unity for central (top) and for 
peripheral (bottom) events, which would imply a scaling close to binary collisions. In central Au+Au 
collisions the protons are significantly lower than the strange baryons. 

Measurements of -Rcp have shown that the Cronin effect cannot account for the relative difference 
between R^a and Rqp between light quark and strange quark particles and among the strange particles 
themselves. It seems that the Cronin effect scales with the hadronic mass rather than the flavor content. 
Thus the Cronin effect for lambdas and protons is almost identical |236l I237j . where the Raa values are 
very different above pt = 2 GeV/c. 

The momentum dependent strangeness enhancement is very similar in central and peripheral col- 
lisions, while the protons are consistent with Raa = 1 in central collisions. The extension of the 
measurement to higher momenta might show if the protons and strange baryons show the same pt 
suppression. It is notable that the A and S Raa values are in agreement with each other, which is 
not the case for the strangeness enhancement factors shown in the previous section. The Q might help 
further to separate the strangeness enhancement from the high pt suppression. 

7 Elliptic flow V2 of strange particles 

The elliptic flow V2 is the azimuthal momentum space anisotropy of particle emission from non-central 
collisions in the transverse direction with respect to the beam direction. Elliptic flow is described by 
the second harmonic coefficient of an azimuthal Fourier decomposition of the momentum distribution. 

Elliptic flow is an observable which is directly related to the initial spatial anisotropy of the nu- 
clear overlap region in the transverse plane, translated into to the observed momentum distribution of 
identified particles if quarks interact and thermalize in the early stage. It is important to understand 
the initial geometry and how it varies with the collision centrality and system size. To compare the 
strength of f 2 of different centralities one has to correct for the initial eccentricity with respect to the 
reaction-plane of the collisions. A review of this experimental measurements can be found in j238] . 

7.1 V2 scaling as a function of momentum 

The V2 results are not unique to strangeness and thus we will only discuss them briefly in the context 
of mass and quark number scaling. Generally, any f 2 discussion needs to distinguish between the low 
Pt regime, where hydrodynamics holds, and the high pt regime, where quark number scaling indicates 
non-bulk production of the particles and their collective behavior. Up to pt = 2 GeV/c the strange 
meson and baryon fit very well in the hydrodynamic systematics, which assumes that the collective 
strength scales with the final hadron mass. Above 2 GeV/c the strange and non-strange particles do 
not scale with the hadron mass but rather with the number of valence quarks. This has been attributed 
to the fact that the partons in the thermal system prior to hadronization are already flowing collectively 
and then simply recombine into a hadron that flows with the added collectivity of all valence quarks. 
This can also be attributed to hydrodynamics but in the partonic system rather than the hadronic 
system, which is one of the signatures of the formation of the deconflned collective partonic phase. 

Figure 1571 shows the experimental evidence for this effect, and an attempt to address the universality 
of the quark number scaling by analyzing the centrality dependence of the scaling properties. 



49 




p^/riq (GeV/c) 



(m^ - m)/n^ (GeV/c^) 



2.5 
2 



> 



X 

^ 0.5 



1 1 1 1 1 

1 Kg 

- ^°oO - 

- 1" 



- e 


1 1 1 1 1 

1 


1 1 1 1 1 

M 

A" ■ 

/ 


• 




s 

1 1 1 1 1 


" (0 

1 1 1 1 1 


0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 

(m^- m)/n^ (GeV/c^) 




(m^ - m)/n^ (GeV/c^) 



Figure 37: Elliptic flow V2 scaled by the number of quarks (Nq) and participant eccentricity (cpart), 
(t>2 / {Nq^ X Cpart)), of identified particles (particle + anti-particle) versus (a) the scaled Pt/N^ and (b) 
(mt — m)/Nq for three centrality bins. For comparison, ideal hydrodynamic model calculations |245] 
are shown as lines in (b). In (c) the same data from (b) is shown, but scaled by the integrated V2 of 
each particle, instead of epart- In (d) the data from (b) scaled by the integrated V2 of all charged hadrons 
is presented. The inset in (d) expands the low rrit region. The error bars only represent the statistical 
uncertainties. All data are for y/s^ = 200 GeV Au + Au collisions [239]. 
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Here the scaled f 2 is normalized first by the participant eccentricity epart derived from a Monte Carlo 
Glauber calculation. The plots (a) and (b) show the doubly scaled quantities from three centrality bins 
as a function of Pt/N^ and (mt — m)/Nq , respectively, which takes out the hadronic mass dependence. 
It is interesting to note that, in both cases, at a given centrality the elliptic flow of all hadrons scales 
as observed for minimum-bias events. Upon dividing by the eccentricity the stronger collective motion 
in more central collisions becomes apparent, which is qualitatively consistent with ideal hydrodynamic 
model calculations, shown as lines in (b). However, there is no universal scaling with eccentricity in the 
STAR data which is in disagreement with the conclusions reached by PHENIX |242] . 

To further clarify the issue, instead of dividing the measured f 2 by the corresponding eccentricity, 
in (c) the averaged V2 as a function of transverse momentum weighted with the measured spectra is 
shown. This scaling seems to work better, however, different hadrons seem to have different values of 
V2, in particular for the top 10% centrality bin at the higher rrit. 

In (d) the doubly scaled f 2 is shown. This time the integrated values of f 2 are extracted from the 
measurements of unidentified charged hadrons in the corresponding centrality bins. In this case the 
scaling seems to work best. It is interesting to point out that at the most central bin (see inset in (d)) 
the values of V2 become negative at low pt for all hadrons. This is most likely caused by the strong 
radial flow developed in central Au+Au collisions, which has also been observed at SPS energies. 

7.2 Partonic versus hadronic flow 

Further evidence for early partonic collectivity is given by investigating in particular the flow of purely 
strange hadrons i.e. the Q and the 0, both of which are expected to exhibit a significantly reduced 
interaction cross section with the surrounding hadronic medium. Figure [39] shows the f 2 for multi- 
strange baryons in STAR. From these early measurements it could already be deduced that the S 
V2 follows the quark number scaling of the non-strange baryons. Recent preliminary high statistics 
measurements by STAR for vector mesons and multi-strange baryons are shown in Fig. |3S] [241j . It is 
evident that they follow in strength and scaling exactly the same pattern than the light quark hadrons. 
Thus a major portion of the collective motion has to be developed in the early partonic stage. 
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Figure 38: Right: Elliptic flow V2 as a function of pt for vr, p (left) and 0, Q (right) in Au+Au minimum- 
bias collisions at y/s^ = 200 GeV from STAR. Open symbols represent results from PHENIX |242] . 
Lines represent a fit inspired by the number of quark scaling |243] . 
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Figure 39: Elliptic flow f 2 of (a) + S+ and (b) VL^ + 0+ for y^s^ = 200 GeV Au+Au minimum-bias 
collisions. The V2 of Kg and A |246] are shown as open symbols, and the results of the fits |243j are 
shown as dashed lines. Hydrodynamic model calculations |244] are shown as dotted lines for K and A 
and as solid lines for S~ and f2~ masses, from top to bottom, respectively |240] . 
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8 Conclusions 



The wealth of data on strange hadron production collected in the recent years at the CERN-SPS and 
RHIC-BNL has allowed to establish a long list of observations on the properties of hot and dense 
hadronic matter produced in heavy ion reactions. The following list summarizes the main findings as 
discussed in the previous sections: 

• The enhancement of the production of strange hadrons per participating nucleon relative to el- 
ementary p-l-p(Be) collisions is a well established experimental fact. The observed enhancement 
increases with the strangeness content of the hadron. 

• A clear onset of the strangeness enhancement with center-of-mass energy has not yet been ob- 
served. In fact, the enhancement of particles with single strangeness increases toward lower 
energies in contrast to naive expectations. 

• Transport models involving only hadronic degrees of freedom are not able to describe the measured 
enhancement of multi-strange baryons. 

• Statistical models generally provide a reasonable description of strange particles yields at all 
energies. The properties of the fireball arc determined by the chemical freeze-out parameters 
and Tch. These turn out to fall on a single freeze-out curve, if the center-of-mass energy is varied. 

• A pronounced maximum of the relative strangeness production is seen around ^/s^ — 7 — 8 GeV 
and the energy evolution of the inverse slope parameters of kaon pt spectra exhibits a drastic 
change in the same region. Both have been interpreted as an indication for a drastic change of 
the properties of the fireball due to the onset of deconfinement. 

• The scaling properties of the (p meson seem to rule out kaon coalescence as the dominant produc- 
tion mechanism. It might be more likely that from SPS energies upwards its origin is to a large 
extend partonic. 

• The transverse mass spectra of S and ft suggest that their kinetic freeze-out is happening in an 
earlier and hotter phase of the fireball evolution. The radial flow observed in these particles would 
thus be mainly partonic by origin and starts to build up in the SPS energy regime already. This 
is corroborated by the substantial elliptic flow of the Q and measured at RHIC, which, due to 
their small hadronic cross section, can only develop in a partonic phase. 

• The measured system size dependence reveals a more complicated pattern as predicted by statis- 
tical models via the transition from a canonical to a grand-canonical ensemble. The contribution 
of a deconflned phase might also alter the system size dependence in a signiflcant way in addition 
to the volume effects in the statistical description. However, the observed behavior can to a large 
extend be explained in a natural way by the core-corona approach. Alternatively, the assumption 
that the strange quark production, in contrast to the A^part scaling of up and down quarks, scales 
with the number of binary collisions works with the notable exception of the meson. 

• Strange particles provide valuable information on the flavor dependence of intermediate and high 
Pt effects, such as high pt suppression. The observed high baryon-meson ratio indicates that quark 
coalescence is dominating over fragmentation at intermediate pt. The observed quark number 
scaling of the elliptic flow at intermediate pt points to a partonic origin of the early pressure in 
the flreball. 
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Starting from the original suggestion of using strange particles as a probe for quark-gluon plasma 
formation, strangeness has turned out to be a versatile tool to study the properties of the matter 
produced in high energy nuclear reactions and has revealed many unexpected effects. Even though the 
interpretation of the observed strangeness enhancement might not be as straight forward as originally 
suggested, it seems there is no consistent explanation without invoking partonic degrees of freedom. 
Models that include only hadronic degrees of freedom are not able to describe the enhancement of 
multi-strange particles. Instead one finds that statistical approaches assuming full equilibration result 
in a much better description. Also the scaling properties of the meson yields and spectra indicate 
that from SPS energies onwards its dominant production mechanism is of partonic origin, i.e. quark 
coalescence, while kaon coalescence cannot explain the observed features. This point of view is further 
strengthened by the phenomenology of radial and elliptic flow of strange particles and follows from the 
fact that particles with low hadronic cross section (0, S, Q) exhibit significant radial flow and that the 
elliptic flow of strange particles scales with the number of constituent quarks. 

Important questions that remain are for instance: Is there a clear onset of strangeness enhancement 
observable when increasing the center-of-mass energy? There are some features in the existing data 
(the maximum of the K^/vr"*" ratio and a sudden change in the energy dependence of the inverse slopes 
parameters of kaon pt spectra) that can be interpreted in this way, but additional support for this 
findings is desirable. Can a unique connection of the enhancement of multi-strange particles to partonic 
equilibration be established or is this a result of other properties of the hadronization processes? What 
is the excitation function of flow phenomena (radial and elliptic) of strange particles and could this 
help to uniquely define the onset of deconfinement? Upcoming measurements at high energies (LHC), 
but also at lower energies (FAIR, RHIC, SPS) will provide important contributions to this strange and 
exciting picture. 
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9 Data compilation 



9.1 Energy dependence 

Figures HQ] and ST] and the Tables [314TT] summarize the current data available on strange particle pro- 
duction in central Pb+Pb and Au+Au collisions. 
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Figure 40: The rapidity densities dN/dy around midrapidity for K"*", K~, Kg, and measured in 
central nucleus-nucleus collisions as a function of y/s^. The systematic errors are represented by the 
gray boxes. 
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Table 3: The rapidity densities dN/dy around midrapidity and the total multiplicities (A^) of K"*" as 
measured in central nucleus-nucleus collisions at different y/s^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((iVpart))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 
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U-0 70 




/ .oOdzU.zidzi.io 
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rjODD/rjyi ( 


|163| 






Au- 


hAu 


n no/ 
U-O/o 




1 1 c:c:_i_n 0/1-1-1 7Q 




Th^QRR /TTOI 7 

rjoDD/rjyi ( 








Au- 


hAu 


n ko/ 
U-0 /o 


o04 


1 1 70-1-0 Q7-I-1 77 
ii. /y=tu.o / zti. ( ( 


OA O-l-H 0-i-Q R 


li/oUz 


1 Rn 
[iou 






Au- 


hAu 


U-4 /o 


oDo 




ZO. (Zti.DltZ.O 


H/OUZ 


[IDl 




R 97 


Pb^ 


-Pb 


U- ( /o 


o4y 


i D . It U . D It U . 


/in 7-i-n 7-1-9 9 

4U. ( ItU. ( ZtZ.Z 


1\T A /IQ 
IN i\4y 


11 9f;i 




7 R9 


Pb^ 


-Pb 


n 7'% 

U- 1 /o 


o4y 


oi O-un 

Zl.ZztU.O_Q g 


oz.yztu.y_3 


1\T A /IQ 
rs i\4y 


h "JAI 
[i^bj 




O. ( O 


Pb^ 


-Pb 


0-7% 


34Q 


on 1 -LA o_Li n 

V7 . J 1— Vy . O— 1 L . U 




NA49 


112^1 




9.2 


Au^ 


hAu 


0-10% 


317 


23.0±4.4(suni) 




STAR 


1247, 




12.3 


Pb^ 


-Pb 


0-7% 


349 


24.6+0.2+1.2 


76.9+2.0+4.0 


NA49 


1125J 




17.3 


Pb^ 


-Pb 


0-5% 


362 


29.6+0.3+1.5 


103.0+5.0+5.0 


NA49 






17.3 


Pb+Pb 


0-4% 




37.1+5.4(suni) 




NA44 


[248 


a) 


17.3 


Pb+Pb 


0-3.7% 




27.5+0.2+1.6 




NA44 


|175| 


b) 


62.4 


Au+Au 


0-10% 




33.6+0.8+2.7 




BRAHMS 


[174; 


c) 


62.4 


Au+Au 


0-5% 


347 


37.6+2. 7(suni) 




STAR 


[122J 




130.0 


Au- 


hAu 


0-6% 


344 


46.3+3.0(suni) 




STAR 


tl22j 




130.0 


Au- 


hAu 


0-6% 




46.2+0.6+6.0 




STAR 


m 




130.0 


Au^ 


hAu 


0-5% 


348 


46.7+1.5+7.0 




PHENIX 


m 




200.0 


Au- 


hAu 


0-5% 


351 


51.3+6.5(suni) 




STAR 


1122] 




200.0 


Au- 


hAu 


0-5% 


351 


48.9+6.3 




PHENIX 


[128 




200.0 


Au- 


hAu 


0-10% 


328 


45.0+0.67+6.75 




BRAHMS 


[129 




200.0 


Au^ 


hAu 


0-5% 






286.0+5.0+23.0 


BRAHMS 


[166) 


d) 
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Table 4: The rapidity densities dN/dy around midrapidity and the total multiplicities (A^) of K~ as 
measured in central nucleus-nucleus collisions at different ^/s^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((A^part))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 



V'^NN 


System 


Cent. 


(iVpart) 


AN/&y 


(AT) 


Exp. 


Ref. 


Comnt. 


(GeV) 


















3.32 


Au+Au 


0-5% 




0.19+0.01+0.03 


— 


E866/E917 


|163| 




3.83 


Au+Au 


0-5% 




0.61+0.02+0.09 


— 


E866/E917 


|163| 




4.29 


Au+Au 


0-5% 




1.26+0.04+0.19 


— 


E866/E917 


[163J 




4.87 


Au+Au 


0-5% 




2.21+0.03+0.33 


— 


E866/E917 


1163] 




4.87 


Au+Au 


0-5% 


354 


2.24+0.05+0.34 


4.14+0.09+0.62 


E802 


[160J 




4.96 


Au+Au 


0-4% 


363 


— 


3.76+0.28+0.38 


E802 


|161| 




6.27 


Pb+Pb 


0-7% 


349 


5.58+0.07+0.11 


10.3+0.1+0.2 


NA49 


1126] 




7.62 


Pb+Pb 


0-7% 


349 


7.8+0.1+0.2 


16.0+0.2+0.4 


NA49 






8.73 


Pb+Pb 


0-7% 


349 


7.58+0.12+0.4 


19.2+0.5+1.0 


NA49 


m 




9.2 


Au+Au 


0-10% 


317 


8.7±2.0(sum) 


— 


STAR 






12.3 


Pb+Pb 


0-7% 


349 


11.7+0.1+0.6 


32.4+0.6+1.6 


NA49 






17.3 


Pb+Pb 


0-5% 


362 


16.8+0.2+0.8 


51.9+1.9+3.0 


NA49 


[125J 




17.3 


Pb+Pb 


0-4% 




21.5±7.5(sum) 


— 


NA44 


[248J 


a) 


17.3 


Pb+Pb 


0-3.7% 




15.4+0.5+1.0 


— 


NA44 


m 


b) 


62.4 


Au+Au 


0-10% 




28.6+0.8+2.3 


— 


BRAHMS 


[174J 


c) 


62.4 


Au+Au 


0-5% 


347 


32.4±2.3(sum) 




STAR 


[122J 




130.0 


Au+Au 


0-6% 


344 


42.7±2.8(sum) 




STAR 


[122J 




130.0 


Au+Au 


0-6% 




41.9+0.6+5.4 




STAR 


|130| 




130.0 


Au+Au 


0-5% 


348 


40.5+2.3+6.1 




PHENIX 


|171| 




200.0 


Au+Au 


0-5% 


351 


49.5±6.2(sum) 




STAR 


|122| 




200.0 


Au+Au 


0-5% 


351 


45.7+5.2 




PHENIX 


1128J 




200.0 


Au+Au 


0-10% 


328 


40.9+0.63+6.13 




BRAHMS 


1129J 




200.0 


Au+Au 


0-5% 






242.0+4.0+19.0 


BRAHMS 




d) 



Table 5: The rapidity densities dN/dy around midrapidity and the total multiplicities (A^) of Kg as 
measured in central nucleus-nucleus collisions at different ^/s^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((A^part))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 





System 


Cent. 


(iVpart) 


AN/dy 


{N) Exp. 


Ref. Comnt. 


(GeV) 














17.3 


Pb+Pb 


0-4.5% 




26.0+1.7+2.6 


— NA57 


|:;!04j 


17.3 


Pb+Pb 




351 


21.9+2.4 


— WA97 


|1111I250| 


130.0 


Au+Au 


0-6% 




33.9+1.1+5.1 


— STAR 


[no] 
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Table 6: The rapidity densities dN/dy around midrapidity and the total multiplicities (A^) of as 
measured in central nucleus-nucleus collisions at different y/s^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((iVpart))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 





System 


Cent. 


(^part) 


dN/dy 


(N) 


Exp. 


Ref. 


Comnt. 


(GeV) 


















4.87 


Au+Au 


0-5% 


340 


0.362+0.065+0.054 




E917 


^182 




6.27 


Pb+Pb 


0-7% 


349 


1.17+0.23+0.38 


1.89+0.31+0.22 


NA49 


|149| 




7.62 


Pb+Pb 


0-7% 


349 


0.94+0.13+0.30 


1.84+0.22+0.29 


NA49 


|149| 




8.73 


Pb+Pb 


0-7% 


349 


1.16+0.16+0.14 


2.55+0.17+0.19 


NA49 


[149] 




12.3 


Pb+Pb 


0-7% 


349 


1.52+0.11+0.22 


4.04+0.19+0.31 


NA49 


[149] 




17.3 


Pb+Pb 


0-5% 


362 


2.44+0.10+0.08 


8.46+0.38+0.33 


NA49 


[149] 




17.3 


Pb+Pb 


0-4% 




2.39+0.21(sum) 


7.6+1. l(sum) 


NA49 






17.3 


Pb+Au 


0-7% 




2.05+0.14+0.25 




NA45 


[181] 


e) 


17.3 


Pb+Au 


0-7% 




2.04+0.49+0.32 




NA45 


[181] 


f) 


62.4 


Au+Au 


0-20% 




3.52+0.08+0.45 




STAR 


[180] 




130.0 


Au+Au 


0-11% 




5.73+0.37+0.69 




STAR 


[183] 




200.0 


Au+Au 


0-5% 




7.70+0.30+0.85 




STAR 


[184] 




200.0 


Au+Au 


0-5% 




7.95+0.11+0.73 




STAR 


[180] 




200.0 


Au+Au 


0-10% 


325 


3.94+0.60+0.62 




PHENIX 


[1851 
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Table 7: The rapidity densities dN/dy around midrapidity and the total multiplicities (A^) of A as 
measured in central nucleus-nucleus collisions at different ^Js^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((A^part))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 





System 




1 AT A 


A AT 1 Aoi 




HjXp. 




Comnt. 




















2.68 


Au+Au 


0-5% 




0.76+0.03 


0.58+0.04 


E895 


[249j 


g) 


3.32 


Au+Au 


0-5% 




3.6+0.1 


5.5+0.3 


E895 


|249| 


g) 


3.83 


Au+Au 


0-5% 




6.8+0.2 


11.6+0.9 


E895 


|249j 


g) 


3.83 


A i A 

Au+Au 




342 




12.0+1.0 


E895 


[131] 




4.29 


Au+Au 


0-570 




10.25+0.3 


16.0+1.0 


E895 


|249] 


g) 


4.87 


Au+Au 


0-5% 




10.0+1.0 




E891 


[143J 




4.87 


A 1 A 

Au+Au 


0-5% 






16.7+0.5+1.7 


E896 


[141] 




8.73 


Pb+Pb 


0-5% 




21.1+0.8+2.1 




NA57 


m 




17.3 


Pb+Pb 


0-5% 




18.5+1.1+1.9 




NA57 


[139J 




17.3 


Pb+Pb 


0-4.5% 


349 


18.5+1.1+1.8 




NA57 


1112J 




17.3 


Pb+Pb 




351 


13.7+0.9 




WA97 


|lllll25()j 




6.27 


Pb+Pb 


0-7% 


349 


13.4+0.1+1.1 


27.1+0.2+2.2 


NA49 


EH] 




7.62 


Pb+Pb 


0-7% 


349 


14.7+0.2+1.2 


36.9+0.3+3.3 


NA49 


[SH] 




8.73 


Pb+Pb 


0-7% 


349 


14.6+0.2+1.2 


43.1+0.4+4.3 


NA49 


[SH] 




8.73 


Pb+Pb 


0-7% 


349 


15.3+0.6+1.0 


45.6+1.9+3.4 


NA49 


[251J 


h) 


12.3 


Pb+Pb 


0-7% 


349 


12.9+0.2+1.0 


50.1+0.6+5.5 


NA49 


[SH] 




12.3 


Pb+Pb 


0-7% 


349 


13.5+0.7+1.0 


47.4+2.8+3.5 


NA49 


[251] 


h) 


17.3 


Pb+Pb 


0-10% 


335 


9.5+0.1+1.0 


44.9+0.6+8.0 


NA49 


[SH] 




17.3 


Pb+Pb 


0-10% 


335 


10.9+1.0+1.3 


44.1+3.2+5.0 


NA49 


[251] 


h) 


130.0 


Au+Au 


0-5% 




17.3+1.8+2.8 




PHENIX 


m 




130.0 


Au+Au 


0-5% 




17.0+0.4+1.7 




STAR 


m 


c) 


200.0 


Au+Au 


0-5% 


352 


16.7+0.2+1.1 




STAR 


[133] 


c) 
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Table 8: The rapidity densities dN/dy around midrapidity and the total multiplicities (A^) of A as 
measured in central nucleus-nucleus collisions at different y/s^^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((A^part))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 





System 


Cent. 


(^part) 


diV/dy 


{N) 


Exp. 


Ref. 


Comnt. 


(GeV) 


















4.87 


Au+Au 


0-12% 




n niQ+o-004 +0.003 

U.Uiy_Q _o.002 




E917 


m 




8.73 


Pb+Pb 


0-5% 




0.44+0.03+0.04 




NA57 


[139] 




17.3 


Pb+Pb 


0-5% 




2.47+0.14+0.25 




NA57 


[I39j 




17.3 


Pb+Pb 


0-4.5% 


349 


2.44+0.14+0.24 




NA57 


111^ 




17.3 


Pb+Pb 




351 


1.8+0.2 




WA97 






6.27 


Pb+Pb 


0-7% 


349 


0.10+0.02+0.01 


0.16+0.02+0.03 


NA49 


m 




7.62 


Pb+Pb 


0-7% 


349 


0.21+0.02+0.02 


0.39+0.02+0.04 


NA49 


m 




8.73 


Pb+Pb 


0-7% 


349 


0.33+0.02+0.03 


0.68+0.03+0.07 


NA49 


[58] 




8.73 


Pb+Pb 


0-7% 


349 


0.42+0.04+0.04 


0.74+0.04+0.06 


NA49 


[251] 


h) 


12.3 


Pb+Pb 


0-7% 


349 


0.82+0.03+0.08 


1.82+0.06+0.19 


NA49 


[sg 




12.3 


Pb+Pb 


0-7% 


349 


1.06+0.08+0.10 


2.26+0.25+0.20 


NA49 


[251] 


h) 


17.3 


Pb+Pb 


0-10% 


335 


1.24+0.03+0.13 


3.07+0.06+0.31 


NA49 


[58] 




17.3 


Pb+Pb 


0-10% 


335 


1.62+0.16+0.2 


3.87+0.18+0.40 


NA49 


m 


h) 


130.0 


Au+Au 


0-5% 




12.7+1.8+2.0 




PHENIX 


m 




130.0 


Au+Au 


0-5% 




12.3+0.3+1.2 




STAR 


m 




200.0 


Au+Au 


0-5% 


352 


12.7+0.2+0.9 




STAR 


1133 
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Table 9: The rapidity densities dN/dy around midrapidity and the total multiplicities (A^) of S~ as 
measured in central nucleus-nucleus collisions at different y/s^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((iVpart))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 





System 


Cent. 


/AT" \ 

(A/part) 


&N jay 




hjxp. 


Hei. Comnt. 


( T\ 

(GeV) 


















3.83 


Au- 


hAu 




342 


— 


0.17+0.05 


E895 


[HI] 


8.73 


Pb+Pb 


0-10% 




1.84+0.16+0.18 


■ 


NA57 


[132] 


17.3 


PbH 


-Pb 


0-10% 




1.91+0.05+0.19 





NA57 


|139| 


17.3 


PbH 


^Pb 


0-4.5% 


349 


2.08+0.09+0.21 




NA57 


[112J 


i ( .O 


PbH 


-Pb 




ooi 


i.O±U.i 




WAy ( 


hill lOKAl 
[1111 UoOj 


6.27 


PbH 


-Pb 


0-7% 


349 


0.93+0.13+0.10 


1.50+0.13+0.17 


NA49 


m 


7.62 


PbH 


-Pb 


0-7% 


349 


1.17+0.13+0.13 


2.42+0.19+0.29 


NA49 


[58] 


8.73 


PbH 


-Pb 


0-7% 


349 


1.15+0.11+0.13 


2.96+0.20+0.36 


NA49 


m 


12.3 


PbH 


-Pb 


0-7% 


349 


1.22+0.14+0.13 


3.80+0.26+0.61 


NA49 


m 


17.3 


Pb+Pb 


0-10% 


335 


1.44+0.10+0.15 


4.04+0.16+0.57 


NA49 


m 


17.3 


Pb+Pb 


0-10% 


335 


1.49+0.08+0.22 


4.12+0.20+0.62 


NA49 


[189] 


130.0 


Au- 


l-Au 


0-10% 




2.00+0.14+0.20 




STAR 


|132) c) 


200.0 


Au- 


(-Au 


0-5% 


352 


2.17+0.06+0.19 




STAR 


[133] 



Table 10: The rapidity densities dN/dy around midrapidity and the total multiplicities (iV) of S"*" as 
measured in central nucleus-nucleus collisions at different ^/s^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((A^part))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 





System 


Cent. 


(A^part) 


AN/dy 


{N) 


Exp. 


Rcf. 


Comnt. 


(GeV) 




















8.73 


Pb+Pb 


0-10% 




0.068+0.021+0.007 




NA57 


[139] 




17.3 


PbH 


-Pb 


0-10% 




0.422+0.023+0.042 




NA57 


[139] 




17.3 


PbH 


-Pb 


0-4.5% 


349 


0.51+0.04+0.05 




NA57 


[112] 




17.3 


PbH 


-Pb 




351 


0.37+0.06 




WA97 


[11111^50] 




7.62 


PbH 


-Pb 


0-7% 


349 


0.05+0.01+0.01 


0.12+0.02+0.03 


NA49 


[sg 




8.73 


PbH 


-Pb 


0-7% 


349 


0.07+0.01+0.01 


0.13+0.01+0.02 


NA49 


[sg 




12.3 


PbH 


-Pb 


0-7% 


349 


0.21+0.03+0.02 


0.58+0.06+0.13 


NA49 


[SH] 




17.3 


PbH 


-Pb 


0-10% 


335 


0.31+0.03+0.03 


0.66+0.04+0.08 


NA49 


[58] 




17.3 


PbH 


-Pb 


0-10% 


335 


0.33+0.04+0.05 


0.77+0.04+0.12 


NA49 


[189] 




130.0 


Au- 


l-Au 


0-10% 




1.70+0.12+0.17 




STAR 


m 


c) 


200.0 


Au- 


(-Au 


0-5% 


352 


1.83+0.05+0.20 




STAR 


[1331 
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Table 11: The rapidity densities dN/ dy around midrapidity and the total multiplicities (A^) of i^^+fi^ as 
measured in central nucleus-nucleus collisions at different y/s^. Also listed are the centrality selection 
(Cent.) and the averaged number of participants ((A^part))- In most cases the statistical (first) and 
systematic (second) errors are given separately, while for a few measurements the quadratic sum of 
both is quoted (denoted by 'sum'). 





System 


Cent. 


(A^part) 


Particle 


dN/dy 


(AT) 


Exp. 


Ref. Comnt. 


(GeV) 






















8.73 


PbH 


-Pb 


0-11% 




n- 




0.085±0.046±0.009 




NA57 


|139) 


8.73 


PbH 


-Pb 


0-11% 




n+ 




0.035±0.020±0.004 




NA57 


|139| 


17.3 


PbH 


-Pb 


0-11% 




n- 




0.259±0.037±0.026 




NA57 


|139| 


17.3 


PbH 


-Pb 


0-11% 




n+ 




0.129±0.022±0.013 




NA57 


ll39j 


17.3 


PbH 


-Pb 


0-4.5% 


349 


n~ 




0.31±0.07±0.05 




NA57 


11121 


17.3 


PbH 


-Pb 


0-4.5% 


349 


n+ 




0.16±0.04±0.02 




NA57 


11121 


17.3 


PbH 


-Pb 




351 


n-- 




0.41±0.08 




WA97 


|lllll^50j 


8.73 


PbH 


-Pb 


0-7% 


349 






0.10±0.02±0.02 


0.14±0.03±0.04 


NA49 


m 


17.3 


PbH 


-Pb 


0-23.5% 


262 


n- 




0.14±0.03±0.01 


0.43±0.09±0.03 


NA49 


m 


17.3 


PbH 


hPb 


0-23.5% 


262 


n+ 




0.07±0.02±0.01 


0.19±0.04±0.02 


NA49 


m 


130.0 


Au- 


l-Au 


0-10% 




n-- 




0.55±0.11±0.06 




STAR 


[135] c) 


200.0 


Au- 


l-Au 


0-5% 


352 


n-- 




0.53±0.04±0.04 




STAR 


11331 
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Comments 

a) Measured in the rapidity range 2.7 < y < 2.9. 

b) Measured in the rapidity range 2.4 < ?/ < 3.5. 

c) The (\.N / dy value is based on the pt extrapolation using the Boltzmann function. 

d) The value of the total multiplicity is based on the extrapolation using the single Gaussian. 

e) Measured in the rapidity range 2.0 < y < 2.4. 

f) Measured via the decay channel — t- e+ + e~ in the rapidity range 2.1 < y < 2.65. 

g) Preliminary data only, shown on Quark Matter 2001 conference [249] . 

h) Data in |251] are not corrected for feed-down from weak decays and are therefore superseded by the 

values given in |58] . 
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List of Variables 



A : Atomic weight 

b : Impact parameter 

c : Velocity of light 

/3 : v/c 

Ps : Surface velocity 

f3± : Transverse expansion velocity 

CT : Lifetime 

dE/dx : Specific energy loss 

dN/dy : Rapidity density 

Es : Strangeness enhancement factor 

(E) I (N) : Energy density / particle density 

7 : Lorentz factor 

7s : Strange quark fugacity 

7q : Fugacity for quark flavor q 

Cpart '■ Participant eccentricity 

/ : Isospin 

J : Spin 

K : String tension 

As : Wroblewski factor 

m : Mass 

me : Electron mass 

mp : Proton mass 

nis : Strange quark mass 

rrit : Transverse mass 

/^B '■ Baryo-chemical potential 

N : Multiphcity 

A^bin : Number of binary collisions 

Alight '■ Number of light {u, d) quarks in a hadron 

-^part) -^p ■ Number of participants 

Aq : Number of quarks of type g in a hadron 

Apc : Number of corona participants 

As '■ Number of strange quarks in a hadron 

n-Q : Baryon density 

Us : Number density of (anti-) strange quarks 

P : Parity 

Pt, P± '■ Transverse momentum 

p : Total momentum 

Q : Q- value of a given reaction 
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Anti-baryon-baryon ratio 
Nuclear suppression factor 
Nuclear suppression factor 
Fireball radius 

RMS width of rapidity distribution 
Radius parameter 
Boost angle p = tanh""^ 
Entropy density 
Strangeness 

Width of rapidity distribution of particle type i (Gauss) 
Inelastic nucleon-nucleon cross section 
Ccnter-of-mass energy of a binary collision 
Center-or-mass energy in nucleon-nucleon system 
Temperature 

Chemical freeze-out temperature 
Thermal freeze-out temperature 
Critical temperature 
Maximum kinetic energy transfer 
Inverse slope parameter 
Equilibration time 

Decoupling times in hydro freeze-out 
Volume 

Surface velocity 
Transverse expansion velocity 
Elliptic flow coefficient 
Yield 

Yield in core 
Yield in corona 
Rapidity 
Beam rapidity 
Atomic charge 
Particle charge 
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